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SUMMARY 


An experimental and analytical study was conducted to assess the importance 
of trailing edge noise as a helicopter rotor broadband noise source. As part of 
this objective, it was desired to develop a generalized noise prediction in which 
the helicopter operating conditions, rotor geometry, and rotor boundary layer de- 
tails were specified. 

Since existing full-scale rotor data are contaminated with many operative 
helicopter noise mechanisms, assessing the importance of trailing edge noise re- 
quired conducting an experiment in which the mechanism could be isolated. This 
resulted in a two phase study in which the first objective was to investigate the 
noise from an isolated two-dimensional airfoil in an acoustic wind tunnel with a 
low turbulence level. The experiment was conducted at close to full-scale heli- 
copter rotor Reynolds numbers using a main rotor blade section to ensure that the 
airfoil turbulent boundary layer characteristics simulated realistic operating 
conditions. Both boundary layer data and acoustic data were obtained for use 
in developing a scaling law for trailing edge noise predictions. The experiment 
also provided a test case for a first principles trailing edge noise theory de- 
veloped as part of the study. 

Results obtained from the isolated airfoil study were extended to the rotat- 
ing frame coordinate system to develop a helicopter rotor trailing edge noise 
prediction. The prediction procedure was tested against selected full-scale heli- 
copter data. Comparisons of the calculated sound pressure levels with the ex- 
perimental data provided as assessment of the trailing edge noise contribution to 
the total helicopter noise spectrum. 

Conclusions resulting from this study can be grouped into two categories. 

In the case of full-scale rotor noise, it was demonstrated that trailing edge 
noise contributes significantly to the total helicopter noise spectrum at high 
frequencies. This noise mechanism is expected to control the minimum rotor noise 
level. A knowledge of the helicopter altitude, speed, angular position, rotor 
tip speed, blade number, rotor chord and span, and local boundary layer thickness 
are sufficient to predict the rotor trailing edge noise contribution. 

In the case of noise radiation from a local blade segment the acoustic direc- 
tivity pattern is predicted by the first principles trailing edge noise theory. 
Noise dependence on local section Idach number varies as Mach number to the fifth 
power in agreement with the scaling law and the fundamental theory. Finally, acous- 
tic spectra from a local blade segment can be predicted using a scaling law which 
includes the Mach number, boundary layer thickness, and observer position. Spec- 
trum shape and sound pressure level are also predicted by the first principles 
theory but the analysis does not predict the Strouhal value identifying the spec- 
trum peak. Accurate. predictions require using measured surface pressure data 
obtained near the airfoil trailing edge. 



INTRODUCTION 


Background 

The noise level at a given angle relative to a helicopter and within a 
specified frequency band represents the sum of contributions of a number of 
separate noise mechanisms. The relative importance of these mechanisms 
depends upon the rotor geometry, aerodynamic characteristics (e.g., tip speed, 
blade pitch), and operating mode of the helicopter (hover, forward flight, 
vertical ascent or descent) . Some of the mechanisms apply primarily to the 
main rotor, some apply to both the main rotor and the tail rotor, and in the 
more complex situation, interactions occur linking the wakes of one rotor 
blade with noise generation by another blade. The simultaneous occurrence 
of the different mechanisms has always complicated both data interpretation 
and attempts to assess the accuracy of theory developed to treat each mecha- 
nism. 


The list of possible rotor noise mechanisms is lengthy and is often 
debated among specialists in helicopter noise research. George (ref. 1) 
provided a comprehensive list of helicopter noise mechanisms. The list can 
be subdivided into periodic and random noise sources both of which are impor- 
tant over the range of helicopter operating conditions. In the absence of 
impulsive noise (ref. 2), George notes that random (or broadband) noise can 
be the primary noise source. For the case of hover, Munch et al. (ref. 3) 
demonstrated using a 2 Hz bandwidth frequency analysis, that broadband noise 
was the prime contributor to PNL for five different rotor systems over a wide 
range of thrust parameters (thrust coefficient divided by solidity) . Study 
of helicopter rotor broadband noise is, therefore, relevant. 

One of the sources of helicopter rotor broadband noise is the rotor tur- 
bulent boundary layer, in particular, its interaction with the blade trailing 
edge. Although there is little dispute that this is an operative mechanism 
in the case of airframe noise, questions exist as to its spectrum shape and 
absolute level on helicopter main and tail rotors. Unfortunately, full-scale 
helicopter rotor data contains other broadband noise sources making it diffi- 
cult to assess the importance of trailing edge noise. The additional mecha- 
nisms include incident turbulence, stalled airfoil noise, and engine noise. 
The objective of the present study was, therefore, to isolate the trailing 
edge noise mechanism and obtain a quantitative understanding of this noise 
source as a function of the relevant parameters. 
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Previous Investigations 


Experimental Studies - Direct experimental studies of helicopter rotor 
trailing edge noise are presently not available. Previous research on this 
noise mechanism was limited to low Mach number two-dimensional isolated air- 
foil tests, full-scale airframe tests, and blown flap noise studies. Although 
these studies did not simulate helicopter rotor operating conditions, they are 
applicable to understanding the trailing edge noise mechanism. 

Previous attempts, such as that of Paterson et al. (ref. 4) to measure 
trailing edge noise from a NACA 0012 two-dimensional isolated airfoil were 
unsuccessful. The airfoil, which was tested in the UTRC open jet acoustic 
wind tunnel, could not be detected above the facility background noise. To 
overcome the background noise sources Schlinker (ref. 5) employed a direc- 
tional microphone system to detect the airfoil. Measurements demonstrated 
that the trailing edge was the dominant source region for airfoil self noise 
generation. Spectrum peak amplitudes were found to vary with approximately 
a fifth power velocity scaling. Comparisons of the measured spectra showed 
reasonable agreement with an existing scaling law prediction developed by 
Fink (refs. 6 and 7) for airframe trailing edge noise. 

Heller and Dobrzynski (ref. 8) circumvented the background noise problem 
by using an aerodynamically clean glider to experimentally investigate 
trailing edge noise. Measured far field acoustic spectra peaked near a 
Strouhal number of 0.02 with the length scale based boundary layer displace- 
ment thickness. Predictions of the radiated noise from the measured surface 
pressure characteristics were found to be satisfactory. 

Yu and Joshi (ref. 9) reported the results of a study using space-time 
correlations to extract isolated airfoil turbulent boundary layer noise 
spectra from facility background noise. The coherent part of the broadband 
noise generation was found to be localized at the airfoil trailing edge in 
agreement with the measurements of Schlinker. Correlation studies concluded 
that the passage of large scale eddies in the boundary layer over the airfoil 
trailing edge generated the noise. This conclusion was important to selecting 
the characteristic length scale needed in the acoustic scaling law developed 
in the present study. A further discussion is given in the section titled 
Airfoil Boundary Layer Characteristics. 

The investigation by Brooks and Hodgson (ref. 10) used a cross-spectrum 
approach to extract the trailing edge noise for an isolated airfoil. Sound 
pressure level dependence on directivity and Mach number were in close agree- 
ment with the theoretical analysis of Howe (ref. 11). 
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Trailing edge noise data has also been obtained from experiments in 
which the turbulent flow over a trailing edge is provided by a wall jet 
(refs. 12, 13, and 14) blowing over the edge of a flat plate. In addition, 
numerous airframe noise studies have provided trailing edge noise data. 

These programs, defined as ”clean airframe” noise studies, provided fly-over 
noise measurements for clean wing and tail surfaces. A summary of the 
various investigations and the resulting conclusions is available in the pub- 
lications of Fink (refs. 6 and 7). The data in reference 13 formed the basis 
of a semi-empirical trailing edge noise prediction procedure which is now 
part of the NASA ANOPP Program. 

Theoretical Studies - Because trailing edge noise is produced by turbu- 
lence (which can be represented by quadrupoles) moving past a sharp trailing 
edge, one analytical approach to the problem is to use the Lighthill acoustic 
analogy (ref. 15). Thus, quadrupole sound sources are placed near the edge 
of a flat plate which produces amplification of the sound if the trailing edge 
is in the near field of the quadrupole. This formulation describes the type 
of model used in the study of Ffowcs Williams and Hall (ref. 16). Subse- 
quently, Crighton and Leppington (ref. 17), Crighton (ref. 18) and Howe 
(refs. 19-21) used this approach in developing trailing edge noise analytical 
models. 

An alternate approach which corresponds to the method employed here, is 
to work directly with the airfoil surface pressures. Since surface dipoles 
are generally more efficient sound radiators than volume quadrupoles, one 
need only consider the quadrupoles to the extent that they induce surface 
dipoles. These dipoles are the source of the amplification mentioned in the 
above paragraph. This approach was used by Chase (refs. 22 and 23) and Chand- 
iramani (Ref • 24) for the case of zero ambient mean Mach number • 

The technique of using the convecting surface pressures was generalized 
by Amiet (refs. 25 and 26) to include an ambient mean flow. Application of the 
Kutta condition produced downstream shed vorticity and also determined the 
magnitude of the induced sound. Brooks and Hodgson (ref. 10) also used the 
convecting surface pressure technique for their comparisons with experiment 
and found good agreement between theory and experiment. If properly treated, 
the first model using Lighthill ’s analogy and the second model using airfoil 
surface pressures should produce identical results, at least to lowest order 
because the surface pressures are determined by the volume quadrupoles. 

Present Study 

Objectives - The overall objective of the present study was to assess the 
importance of trailing edge noise as a helicopter rotor broadband noise source. 
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As part of this objective it was desired to develop a generalized noise 
prediction in which the altitude, speed, angular position of the helicopter, 
rotor tip speed, and rotor boundary layer details were specified. Hover was 
to be included as a special case. 

Since existing full-scale and model scale rotor data are contaminated 
with other operative helicopter noise mechanisms assessing the importance of 
the turbulent boundary layer broadband noise necessitated conducting an experi- 
ment in which this mechanism could be isolated. This resulted in a two-phase 
program in which the first objective was to investigate the noise from an 
isolated two-dimensional airfoil in an acoustic wind tunnel with low incident 
turbulence levels. This simulation technique is shown in figure 1. In addi- 
tion to obtaining acoustic data, the isolated airfoil study was required to 
provide turbulent boundary layer information at the airfoil trailing edge. 

The aerodynamic and acoustic data from the experimental study were then used 
to develop scaling relationships for trailing edge noise radiation from a 
two-dimensional airfoil segment. This effort was aided by the empirical 
formulation of Fink (ref. 6) and a separate analytical formulation developed 
as part of the present study. Both the scaling law and fundamental theoreti- 
cal approach were compared with the several isolated airfoil studies to assess 
the accuracy of each approach. 

The objective of the second phase of the study was to extend the station- 
ary two-dimensional isolated airfoil results to the rotating frame case to 
provide full-scale rotor broadband noise predictions. Based on the generalized 
noise prediction procedure, the importance of trailing edge noise could then 
be assessed. This objective was achieved by comparing rotor trailing edge noise 
predictions with the total noise spectrum measured during a full-scale heli- 
copter flyover when all noise mechanisms are present. 

Approach - The first phase of the study was conducted using an untwisted 
segment of a helicopter main rotor. The isolated airfoil section was mounted 
between parallel sideplates (to provide two-dimensional flow) and tested in an 
acoustic wind tunnel (see fig. 2). Freestream Mach number was varied from 
M = 0.1 to 0.55 providing Reynolds numbers close to full-scale values. Noise 
and mean velocity profiles at the airfoil trailing edge were measured for a 
range of test conditions. A directional microphone was used to separate 
trailing edge noise from other facility background noise sources. Directivity 
information was acquired by locating the directional microphone system at 
different viewing angles relative to the airfoil. Geometric angle of attack 
was varied over a range of a = 0° to 12° to determine the effect of this param- 
eter on the boundary layer development and the resulting noise radiation. 

Test results were analyzed to determine scaling law relationships for 
trailing edge noise radiated in the overhead plane of the two-dimensional 
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isolated airfoil. Included in the scaling law evaluation was the previous 
boundary layer noise data reported by Schlinker (ref. 5). and the data of 
Brooks and Hodgson (ref. 10). These different experimental results provided 
a data base with widely varying freestream Mach number (M = 0.1 to 0.55), 
airfoil chord (c = 21 cm to 61 cm), and airfoil profiles to test the generality 
of the scaling law. 

During the second phase of the study results of the two-dimensional 
isolated airfoil investigation were extended to the rotating blade case to 
provide predictions for the full-scale flyover conditions. The transforma- 
tion from stationary coordinates to rotating coordinates assumed that span- 
wise segments of the rotor blade, shown in figure 1, radiate boundary layer 
trailing edge noise with the same spectrum and directivity as an isolated 
airfoil. The local relative velocity of the rotor blade segment determined 
the boundary layer details and the resulting absolute noise levels and spec- 
trum shape. The generalized boundary layer details employed in the noise 
prediction were determined by flat plate boundary layer theory. 
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LIST OF SYMBOLS 

Transmitted acoustic rays; see figure 14 

Decrease in directional microphone response due to scattering 
by turbulence 

Apparent source to shear layer separation distance in flow 
plane; see Appendix C 

Semichord = c/2 

Apparent source-to-shear layer separation distance in vertical 
plane; see equation (C.12) 

Sound speed 

Airfoil chord 

Directional microphone reflector aperture diameter 

Directivity factor defined by equation (26) 

Normalized directivity factor defined by equation (32b) 

Difference between airfoil surface pressure spectra and flat 
plate surface pressure spectra 

Focal point microphone diameter 

Fresnel integral combination defined by equation (9) 

Function describing scaling law spectrum; see equation (67) 
Spectrum function defined by equation (22b) 

Frequency 

Frequency parameter defining scaling law spectrum peak; see 
equation (67) 

Directional microphone gain; see equation (A.l) 

Normalized airfoil response function; see equation (11) 

Directional sensitivity of reflector 
Source-shear layer separation distance; see figure 13 
Unit vectors in x, y, z directions, respectively 
/-I 

Constant in trailing edge noise theory; see equation (74) 

Constant in scaling law for trailing edge noise; see equation (68) 
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Wavenumbers for x and y dimensions 

Effective lift function defined by equation (13a) 

Spanwise length scale of surface pressure; see equation (19) 

Freestream Mach number 

Spanwise flow Mach number 

Turbulence convection Mach number 

Chordwise flow Mach number 

Flight Mach number in rotor plane 

Mach number of source relative to fluid 

Mach number of flow relative to rotor hub 

Tip Mach number 

Axial flow Mach number 

Overall sound pressure level 

Vector and unit vector from retarded source position to observer 
Magnitude of pressure wavenumber component 
Surface pressure jump induced by edge 
Sum of and p^ 

Surface pressure incident on edge 

Mean-square of focal point microphone sound pressure level 
Mean-square of omnidirectional microphone sound pressure level 
Dynamic pressure of freestream 

Acoustic source to reflector; distance for on-axis source 
Reflector radius of curvature 

Acoustic source to reflector distance for off-axis source 
Reynolds number 

Critical Reynolds number for transition to turbulence 
Revolutions per minute 

Fourier transform of surface pressure wavenumber spectrum; see 
equation (16) 

Fourier transform of R^q; see equation (17) 

Source to observer distance 

Radius of microphone array in retarded source coordinate system 
Focal point microphone distance from reflector surface 
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Designates apparent source in horizontal plane; see figure 47 
Strouhal number 

Designates apparent source in vertical plane; see figure 47 
1/3 octave band sound pressure level 
Sound spectrum 

Instantaneous sound spectrum as a function of azimuthal angle 

Surface pressure spectrum for a flat plant and an airfoil 
Airfoil span 

Distance along ray path between source and shear layer 

Boundary layer trip thickness; also airfoil trailing edge thickness 

Propagation time of sound to reach observer 

Time 

Distance along ray path between apparent source, Sy, and shear layer 
Freestream velocity 

Convection velocity of surface pressure pattern 

Local velocity in boundary layer normalized by freestream velocity 
Diffraction pattern half-width 

Cartesian coordinates (also, x is the distance defining the off- 
axis source position; see figure 12) 

Coordinate system defined by equation (48) 

coordinate system rotated about z^ by y 

Xp coordinate system rotated about y 2 by y 

Observer position in 3c coordinate system 

Source position defined by equation (46) 

Retarded coordinates of rotor hub 

Sideline distance 

Angle of airfoil out of flow plane; also geometric angle of attack 

Azimuthal angle of rotor; also intermittency parameter 
Boundary layer thickness 
Boundary layer displacement thickness 
Small parameter 
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Parameter defined by equation (7) 

Specific value of ? given by equation (15) 

[ (1-M cose^)^ - cos2e^]^/^ 

Diffraction parameter - it Df sin 4 >q/Cq; also, nondimensional 
parameter controlling reflector response to off-axis sources 

Parameter defined by equation (15) 

Observer angle defined in figure 26 

Observer angle for rotor problem defined in figure 27 
Propagation angle of acoustic ray inside open jet (see figure 13) 
Angle defining retarded source position 

Propagation angle of acoustic ray in the absence of flow (see 
figure 13) 

Propagation angle of transmitted acoustic ray outside open jet 
(see figure 13) 

Wavenumber given by equation (5) 

Parameter defined by equation (7) 

Kinematic viscosity of air 
Density of air 

Nondimensional frequency, w = cu6*/U 


THEORETICAL FORMULATION OF THE TRAILING EDGE NOISE MECHANISM 


Noise Model for Isolated Airfoil in Rectilinear Motion 

The foundation for the trailing edge noise model developed in the present 
study is described in references 25 and 26. The model assumes a frozen sur- 
face pressure pattern which convects downstream at a Mach number M^. past the 
trailing edge. Upon approaching the trailing edge the pressure jump across the 
airfoil is forced to zero to satisfy the Kutta condition. This gives rise to an 
induced pressure field which propagates away from the trailing edge region. 

The surface pressure jump, Ap^^due to the convecting turbulence is decomposed 
into spatial Fourier components of the form 


where the time dependence exp (loot) has been removed. Here x is the chordwise 
coordinate, y the spanwize, and z the normal. The airfoil lies in the z = 0 plane 
and X = 0 designates the trailing edge. The parameter A is related to the chord- 
wise gust convection speed by 


X = io/U^ ( 2 ) 

Satisfying the Kutta condition at the trailing edge then gives rise to an 
induced surface pressure jump given by equations (12) and (32) of reference 27 
which is 


Apj (x,y ,2 = 0) = Apf(f,y) 




r+C 


(3) 


Here r (= -x) represents the distance between a point, x, on the airfoil and the 
trailing edge. 


Also 

_ 2 . 2,«2 
5 = /I. - ky />3 

(4) 

/3^ s 1 - 

(6) 


/X s 

(5) 

H tu/u 

(7) 


where M is the freestream Mach number. This equation was derived for a semi- 
infinite flat plate airfoil with a trailing edge but no leading edge. A correc- 
tion for the leading edge could be added to the solution in an iterative scheme 
as described in reference 27, but this will not be needed here. 
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Equation (3) was given in reference 27 for an airfoil surface pressure rather 
than a surface pressure jump. However, since the problem treated in reference 27 
is antisymmetrical abov^ and below the airfoil, the results apply as well to a 
pressure jump. For the present problem writing the equation in terms of a sur- 
face pressure jump is necessary since the problem is not antisymmetrical. Far 
upstream of the trailing edge the surface pressure is produced by volume quadru- 
poles and induced surface dipoles which double the pressure on the side containing 
the turbulent eddy and cancel it on the other side. Downstream of the trailing 
edge the pressure is produced by the frozen quadrupoles resulting in a pressure 
field whose magnitude is one-half of the pressure upstream. 

Reference 26 mistakenly says that the pressure is forced to zero by the 
Kutta condition at the trailing edge, but it is really the pressure jump that is 
being forced to zero. This misstatement has caused some confusion to readers. 

The equations in reference 26 are for the pressure jump, however, and the pre- 
dictions therein are not affected by this misstatement. 

Substituting the assumed incident surface pressure (or surface pressure 
jump) given by equation (1) into equation (3) gives for the induced pressure jump. 


where 


Apj(x,y, 2 = 0) = - Pq e |i - (i + i)E* [r(X+ftM + C)]} 


E(x) s 




_dL 


( 8 ) 


(9) 


is a combination of Fresnel integrals, and the star on E represents the complex 
conjugate. 

The total surface pressure jump is the sum of incident and induced pressure 
jumps. The surface loading is equated to a dipole strength which is used to cal- 
culate the noise generated. In the original derivation given in reference 25, the 
incident pressure jump was erroneously omitted in calculating the noise; reference 
26 gives the corrected results, including the incident pressure jump. However, 
equation (1) is not an ideal representation of the incident pressure, as noted in 
reference 26, since it discontinuously arises at the leading edge of the airfoil, 
rather than growing gradually with the growth of the boundary layer. This would 
give rise to an erroneous leading edge term in the noise calculation. Thus, 
equation (1) will be modified slightly with an exponential convergence factor, 
exp (cXx), so that the incident pressure at the trailing edge (x = 0) is unchanged 
from equation (1), but is small at the leading edge (x = -c) . This will be the 
case if 
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eXc » I 


(10a) 


(10b) 


( 11 ) 


Also, e will drop out of the result if 

€ « I 

Thus, the total surface pressure jump is taken to be 

APs = Poe‘‘‘^’‘"''y^’{e*’^*-l + (i + i)E*[r(X + ^M + C)]} 

s Pog(x, Kx, X, ky) -C < x < 0 

The last identity defines the normalized response function, g. 


For an airfoil in rectilinear motion the far-field sound can be calculated 
in the same manner as described in reference 28 for leading edge noise. This 
was the method used for trailing edge noise described in references 25 and 26. 
The normalization of g in reference 28 differs from that in equation (11) above 
by the factor 2 Also, the factor U in equation (18) of reference 28 

should be for a surface pressure convecting at other than the freestream 
velocity. Taking account of this factor, equation (15) of reference 28 for the 
far-field sound spectrum of an airfoil in rectilinear motion becomes 


2 

Spp(x,w) = ^ 2wa^ ^ |^(x,Kx,X,Ky)| 4>qq(X,K^ 


( 12 ) 


where 


1^1 ^1/ 9(x„.k„X,k 

-c ^ b 

Ky S wy/Cq a 


(13a) 

(13b) 

(13c) 


s is the airfoil span and 4>qq is the wavenumber spectrum of the airfoil surface 
pressure produced by the turbulence. The function ^ can be thought of as an 
"effective lift". For an observer directly above the retarded source position 
(x = 8z, y = 0) , all points on the airfoil will be at equal retarded distances 
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from the observer. For this case the exponential phase factor in equation (13a) 
drops out and 2 becomes the normalized lift. 

Substitution of equation (11) for g into equation (13a) and assuming equation 
(10) holds gives f or 


V Mx/a + Co 


(l + i)E*(2(X+/ZM + Co))] 
(l + l)E*(2(/ix/a+Co))| 


(14) 


where 

®= X + ^(m-x/ct) 


_ = Myx^+/32z2 /o- 
" y ” 


(15) 


and where the bar over X, p and ^ indicates normalization by the semichord b. 


The only remaining unknown in equation (12) is the surface pressure spectrum 
of the turbulence. Since the wavenumber spectrum $qq and the cross-correla- 
tion function Rqq are Fourier transforms of one another 


The quantity Rqq (X,y) can be related to the cross-spectrum that would be 
measured by two probes separated by a distance y. Because the turbulence is 
assumed to be frozen and convects at a speed U^, it can be described by a corre- 
lation function Rqq (x-U^t,y). For fixed t = t^^, the function Rqq (x-U^t,y) 
represents a spatial cross correlation while for fixed x = x^ the function 
represents a timewise cross correlation, i.e., the cross-correlation function 
obtained from two fixed surface probes as the fluid convects by. The Fourier 
transform of the first case with respect to x gives the wavenumber spectrum Rqq 
while the Fourier transform of the second case with respect to time gives the 
frequency spectrum Sqq. The two are related by the factor as noted below 


Rqq(X,y)=^ / Rqq (x-Upt.y) e~'^*dx 
-00 

^ ^ /°°Rqq^Uc(T-t),y) dr 


(17) 


-00 


I 
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^qq (***» y) 



For an observer in the y = 0 plane, Ky = 0 and equations (16) and (17) give 


^^(X.O) = ^ /®S («.y„)dy„ ( 18 ) 

0 

where the fact that Sqq (o),y) is symmetric about y = 0 was used. The spanwise 
correlation function ly (oj) can then be defined as 

(19) 

0 

For y 0, equation (18) would contain a factor exp (-i Ky yo) under the integral, 
and equation (19) could not be used to replace this integral. However, here 
l>qq(Kj^,Ky) will be assumed to be independent of Ky; i.e., the assumption is made 
that 


4*qq(X, Ky) » (X,0) 


( 20 ) 


Then equation (12) becomes 

iy (o») Sqq (cu, 0) (21) 

The assumption given by equation (20) does not eliminate the angle from the direc- 
tivity prediction. It is only a statement of the details of the boundary layer 
surface pressure distribution, and although this distribution affects the direc- 
tivity of the noise generated, it is believed that the other factors such as 
spanwise correlation length, which are retained in the analysis, are more 
important . 

As in ref erences 25 and 26, Sqq (w,o) will be approximated using measurements 
for flat plate boundary layers. Thus, the surface pressure spectrum data of 
Willmarth and Roos (ref. 29) (also presented as fig. 7-39 of ref. 30) for a tur- 
bulent boundary layer flow of density and freestream speed U can be approxi- 
mated by the expression 



Sqq(o;,0) = (-^ ^ 2 x 10*® F(CJ) 

(22a) 

F(QJ) s (I + © + 0.217©^ + 0.00562 ©V 

(22b) 



1 F(©)d©«3.22 

(22c) 

-00 
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where & = 6 and 6 is the turbulent boundary layer displacement thickness 

given approximately by (ref. 31) 


S*/C = 0.047 Re° ^ (23) 

and Re^ is the Reynolds based on chord. The spanwise correlation length can be 
found by graphical integration of the curve for B given in figure 13 of 
reference 32. This gives the value 

iy(cu) « 2.1 Uc/w (2A) 

where is the convection velocity of the turbulence. In general is a 
function of frequency, but will here be set equal to the average value 

Uc = 0.8U (25) 

It is realized that the expressions given by equations (20) to (25) are approx- 
imations which can lead to errors. This is especially true since they are intended 
for infinite flat plate boundary layers whereas for the present case there is a 
trailing edge which can affect the flow in the immediate vicinity. Difference 
between flat plate surface pressure spectra (in the absence of a trailing edge) 
and isolated airfoil surface pressure spectra are discussed in the section titled 
Trailing Edge Noise Predictions. Experimental measurements are used to quantify 
the differences between these configurations and predict the far field acoustic 
spectra for the isolated airfoil and helicopter rotor blade configurations 
investigated in the present study. 

Equations (13) to (25) are the basis for the theoretical prediction method 
presented herein. (A semi-empirical prediction is given later.) It is inter- 
esting to examine certain of these expressions in limiting cases for comparison 
with work done previously. Because Ky is set equal to zero in $pp(K^,Ky). the 
only remaining factors contributing to the directivity function D (0,(()) are 

0(ee»=(^) 

where 0 is the angle between the downstream x axis and the source observer 
line and (J> is the angle between the xy plane and the plane containing the x axis 
and the observer (see figure 3). Thus, 

X = r cos 6 

y = r sin 0 cos <#> = yg (27) 

z = r sin 0 sin <#> = Zg 
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The retarded coordinates tg and 0g are identified by the subscript e and can be 
shown to be (see e.g., reference 39) 



X = fe (M + cos 6 q ) 


o-= re(i + Mcos 6 q ) 

For the case of large chord (or high frequency) 

Xb » I 

K^b » I 


and using the relation 


equation (14) becomes 


fjib » I 


E*(x) ^ 

' x—co 2 


1 ^ 1 -i 


X /xM Co 

^x/cr + Co 


( 28 ) 


(29) 


(30) 


(31) 


In the y = 0 (<|) = 7t/2) plane = p so that from equations (26), (27), (28), 
and (31) 

,-2 

D (©e. ■^/2) — 2 sin^ (0g/2) (l + M COS 6 ^)'' [l + (M- Mg) cos ©gj <32a) 

where 

D{6q, 4>) ^ 0{de,<p) / D (32b) 

A result equivalent to this has been obtained by other investigations. For 
example, Howe (ref. 11) obtains equation (32) but with one additional factor since 
he allows the shed vorticity to convect downstream at a speed unequal to the 
stream speed. If these two speeds are equated, as is assumed here, Howe's result 
is the same as equation (32). Goldstein (ref. 42) also obtains a directivity 
expression for trailing edge noise, but in the presence of a shear layer. If 
equation (3.44) of Goldstein's paper is corrected for shear layer refraction in 
the manner described by Amiet (ref. 34), the result is identical to equation (32) 
above. The details of this comparison with Goldstein's results for both leading 
and trailing edge noise is to be published in J. Sound and Vibration. Additional 
discussions of Goldstein's analysis are given in the section titled "Experimental 
Assessment of Trailing Edge Noise Radiation From a Two-Dimensional Airfoil Section." 

For low frequency or small chord/wavelength 





(33) 
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Thus, for an observer in the y = 0 plane 



sin^^e 

(i + M cos Qof [i + (M - Me) cos 0^ 


(34) 


This gives a compact dipole like behavior in that the directivity varies as 
sin^ 9g rather than sin^ equation (32). Equation (34) is not precisely 

the directivity of a dipole which would be 



sin^ge 

(l+ M cosge)^ 


(35) 


The reason for this difference is that the assumptions given by equation (10), 
which was used in deriving equation (11), no longer hold for very low frequency. 
This could be rectified by retaining in equation (14) the exponential term from 
equation (11). However, this is not considered necessary as trailing edge 
noise is basically a high frequency phenomenon localized near the trailing 
edge. As frequency decreases causing the wavelength to approach the chord 
length it becomes improper to describe the noise as coming from a small region 
near the trailing edge because the leading edge gives a contribution. Thus, at 
low frequency it becomes less meaningful to consider trailing edge noise as a 
discrete mechanism. Verification of the high frequency characteristics of 
trailing edge noise is provided by the good agreement between the directivity 
pattern of equation 32 and the measurements obtained in the present study 
(see subsections titled Directivity Assessment and Comparisons with Previous 
Investigators) . 


Rotating Blade Noise Prediction 

Equation (21), together with the associated equations for £ and Spp, gives 
the far-field trailing edge noise of an airfoil in rectilinear motion. This 
formulation must now be transformed to the rotating frame to provide a propeller 
or helicopter rotor trailing edge noise prediction. The method employed assumes 
that the spectrum of a given blade segment is, at any particular instant, given 
by equation (21) . This "instantaneous spectrum" varies as the rotor moves about 
the azimuth. To find the final spectrum, this "instantaneous spectrum" is 
averaged around the azimuth, together with a weighting factor which accounts for 

retarded time effects. This is the method used in references 35-37, and is a 

standard technique discussed further by Bendat and Plersol (ref. 38). 

The present formulation for trailing edge noise is somewhat simpler in com- 

parison to studies in references 35-37. In the previous studies, the problem of 
rotor- turbulence interaction was modeled analytically, and this involved treating 
blade-blade force correlations since more than one blade could intersect a given 
eddy. For trailing edge noise, the mechanism for producing blade-blade correla- 
tion is absent and the noise from several blades is determined by the noise from 
a single blade multiplied by the blade number. 


18 


The present analysis will, thus, be concerned with calculating the average 
spectrum produced by a single blade segment. Figure 4 shows the geometry of 
the problem. The origin of the x, y, z axis is fixed to the rotor hub with z 
being the rotor axis. The axial component of flow is in the negative z direction. 
The observer is assumed to be in the x~z plane at a distance r from the hub and 
at an angle 6’ with the z axis. The nonaxial component of flow is at an 
angle ij; to the y axis, pointing inward as shown. The observer is fixed relative 
to the rotor hub, and the x, y, z coordinate system. 

The blade segment is assumed to be a flat plate making an angle a with the 
x-y or azimuthal plane where 


M*+MxCOS(y+'^) 

cot a = — — I 

M2 


(36) 


where = azimuthal Mach number of the blade (relative to the fixed observer, 
not the fluid). Equation (36) assumes that the airfoil is aligned with the flow, 
implying zero steady loading. Because of the skewed inflow, there will be a 
spanwise component of flow over the blade where 

Mq= -M jSin (y+t)(tsin y + t cos y) ^37^ 

The chordwise component M^j^of flow Mach number is 


^ch" (y + ’^)] ("t sin y + 1 cos y) + M2 


(38) 


The expression given by equation (21) for the far-field sound of an airfoil 
in rectilinear motion is written in terms of the present airfoil position. In 
order to be able to apply this relation to the rotating airfoil segment, the 
equivalent of the present position must be calculated for the airfoil segment. 

This "present position" is the position of the airfoil segment relative to the 
fluid if it were to move along a rectilinear path after sound generation until 
that sound reached the observer. An alternative method would be to calculate the 
retarded source position and rewrite equation (21) in terms of the retarded source 
position. Both methods should be of equal difficulty and should give the same 
results; the former will be used here since it is the method used in the previous 
references. 

Both methods begin with a calculation of the retarded source position of the 
airfoil segment (or the rotor hub since distances on the order of the acoustic 
source size are not important when calculating the source position for a far- 
field observer). If at time t = 0 a marker is placed in the fluid and a burst 
of sound is emitted, at time t = Tg, when the observer hears the sound burst, the 
marker will move to position Xg. The observer is at 
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Xq = r(t sin 0'+ k cos 6') 


(39) 


which corresponds to a distance tg from the retarded source point Xg. Thus, 


= (r sin 9'- Xg)^ + + (r cos 6'- Zg)^ 


(40) 


Also, 


Tg - fe/Cg 


(41) 


Since and are given by 


Mj = -Mjk 


(42) 


Mf = -Mj(l sin t + j cos 

the retarded source coordinates are 

Xg = -Tg Cg Mf sin t = - Mf Tg sin f 

yg = -Mjrg cos»^ (43) 

Zgs-M^rg 

Substitution of equation (43)into equation (40) leads to the following result for rg 


= " 2 [Mg COS © + yi - Mg^ sin^ © j 


(44) 


where 


Mg s M| + Mj 
Mg^ = Mf^ + m/ 


(45) 


and Mg cos 0 = M^ sin \p cos 6 ' + Mg cos 0*. 

By calculating (Mg/Mg) • (-x^) , 0 is shown to be the angle between the convec- 
tion Mach number Mg and the vector linking the observer and the source, -Xg. 
Substitution of equation (44) into equation (43) determines the retarded source 
position Xg. 
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The "present** source position Xp is found by adding to the retarded source 
vector a vector giving the chordwise airfoil movement during the time T^. The 
spanwise component of fluid velocity is ignored here since it was ignored in 
deriving equation (21). This is possible because in calculating the noise of an 
infinite span airfoil in rectilinear motion, any spanwise flow velocity can be 
eliminated by a coordinate transformation. The only purpose of calculating 
the ’’present** source position is to find the appropriate coordinates for use in 
equation (21) . This position is 

ip " (46) 


Introducing equations (43) and (38) into equation (46) gives 

ip/fe = cos y sin (y+ V')] 

+t [m^ cos y- M| sin y sin (y+ '/')] 


(47) 


For a coordinate system at this ’’present’* source position the observer has 
coordinates xj_ given by 




Xo 




(48) 


Equation (21) is given in a coordinate system in which the flat plate airfoil 
lies in the x-y plane with the span along y. The x^ coordinate system above 
must be rotated to the same orientation with respect to the airfoil. Rotation 
of the xj system about the axis by an angle y gives the X 2 system with y 2 
now pointed along the airfoil span 

Xg = K, sin y- Yj cos 7 


Vg = X| cos y + Yj sin X 


^2 = 


(49) 


By rotating the X 2 system about the y 2 axis by an angle ot gives the x^ system 
with X 3 along the chord pointing from leading to trailing edge and 

X3 = Xg COS a - Zg sin a 

ys=y2 (50) 

Z3 = XgSin a + Zg cos a 

The x^ system is the appropriate system for calculation of the observer position 
needed in equation (21). Introducing equations (47) to (49) and equation (39) into 
equation (50) gives for the observer position 
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where 


and 


X3 = Tg COS a - Tq COS 

Y3 = Xq cos y+ Mf fg Sin(y+'/') 

Z3 = (Xq sin y+fgMf) sin a + Zq cos a 

cos = cos 0' sin a - sin 6'sin a cos a 
^ 0 ^ = Xo^ + Vo^ + li 


( 51 ) 


(52) 


(5-3) 


Here <1^ is the angle between the axis and the actual airfoil position (not the 

"present" position defined above); this can be seen by substituting x^ for x^ in 
equation (49) and proceeding to the calculation of X 3 in equation (51). 

Equations (51) to (53) are the appropriate "instantaneous observer coordinates" 
to use in calculating the sound spectrum given by equation (21). The angle y - Qt 
is a function of time. In order to calculate the time averaged spectrum, an 
average around the azimuth must be determined. This average must account for 
both the Doppler shifting of frequency as the airfoil segment moves relative to 
the observer, and the retarded time effects as the blade moves around the azimuth. 
The latter correction is needed because the blade spends different amounts of 
time, in an acoustic sense, at each azimuthal location. As noted in reference 35, 
the proper azimuthal weighting is the factor (jl)/(a)q where co is the frequency of the 
airfoil forces and is the Doppler shifted frequency. The azimuthally averaged 
spectrum is then 


Spp(x,cj) = 


27T J OJq 


^0’ ^ ) d y 


(54) 


The spectrum Spp (x, o)^, y) for the case of relative motion between source and 
observer is related to Spp (x, o), y ) with no relative motion by the Doppler factor 
0 )/ 0 )^; i.e., as shown in reference 39 . 


Spp(^» ^0*^^ * Oi)o 


(55) 


Thus, equation (54) and (55), together introduce a factor of (w/cjo^) 
The Doppler factor is shown in reference 39 to be 


^ Mt OS 


(56) 


I-M 50 S 
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The final result from the azimuthally averaged spectrum is given by 
equations (54), (55), and (59) with Spp (x, lo, y) in equation (55) given by 
equation (21). The x, y, z coordinates to be used in equation (21) are 
X 3 , y 3 » Z 3 , respectively, given by equation (51). 
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DESCRIPTION OF THE EXPERIMENT 


Acoustic Research Tunnel 

Operating Characteristics - The experimental study was conducted in the 
UTRC Acoustic Research Tunnel. A detailed description of the facility is given 
in reference 40. The tunnel, shown schematically in figure 2, is an open-circuit, 
open-jet design. The inlet is provided with a high length-to-diameter ratio 
honeycomb section and a series of turbulence suppression screens. These features, 
in conjunction with a large tunnel contraction, provide a spatially uniform, 
temporally steady flow with a controlled test section turbulence level of approxi- 
mately 0.2 percent. Turbulence generators and grids can be inserted upstream of 
the nozzle to generate wake profiles and a range of turbulence levels in the test 
section. 

The open jet test section is surrounded by a sealed anechoic chamber 4.9 m 
high, 5.5 m long (axial direction), and 6.7 m wide. Downstream of the test sec- 
tion the airflow enters a diffuser by way of a circular collector that has anechoic 
treatment on its flow impingement surface. The diffuser operates unstalled and 
is thus not a major source of background noise. To avoid tunnel fan noise from 
propagating upstream into the anechoic chamber, a Z-shaped muffling section with 
two right-angle bends and parallel treated baffles is located between the diffuser 
and the fan. The 1100 kW centrifugal fan exhausts to the atmosphere through an 
exhaust tower. 

Tunnel speed is determined from total pressure measurements at the contrac- 
tion inlet and static pressure measurements within the sealed anechoic chamber. 
Since losses are confined to the boundary layer, total pressure upstream and 
downstream of the contraction are predicted and have been verified, to be equal. 

The test section velocity has been shown to be temporally steady. 

Open Jet Geometry and Test Section Arrangement - A rectangular 0.79 m by 
0.53 m test section provided two-dimensional flow conditions for the helicopter 
main rotor blade segment. The blade segment was mounted vertically in the 0.53 m 
direction, between the two horizontal sideplates shown in figure 5. The sidewalls 
provided a uniform spanwise mean flow and eliminated the problem of additional 
noise generation by incident turbulence had the airfoil extended through the 
thick turbulent open-jet shear layer. The airfoil mounting plates on each side- 
wall permitted angle-of-attack variations. 

Triangular tabs previously used to suppress the edge tone phenomenon (ref. 40) 
between the inlet nozzle and collector were not needed. This was because the 
as)nnmetry produced by the rectangular jet flowing into the circular collector 
destroyed the flow field spatial coherence responsible for edge tone generation. 
Removal of the edge tone suppression tabs was investigated in the earlier trailing 
edge noise study of Schlinker (ref. 5). Overall background noise was reduced by 
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approximately 5 dB by removing the tabs. Unfortunately, this improvement was 
not sufficient to detect the airfoil trailing edge noise over the facility 
background noise. As in the present study, a directional microphone was needed 
to measure the airfoil self noise. 


Model Airfoil 

The objective of the experimental phase of this study was to measure tur- 
bulent boundary layer noise at Reynolds numbers close to full scale. For 
this purpose, a 0.53 m spanwise segment of a helicopter main rotor was employed. 
(Note that all acoustic spectra measured in the present study were converted 
to the acoustic radiation from a unit spanwise length (1 ft or 0.3 m) of the 
model airfoil.) The airfoil segment, shown in figures 5 and 6, is part of a 
larger main rotor assembly used for aerodynamic response tests conducted at 
UTRC (ref. 41). The test airfoil with its 0.41 m chord permitted varying the 
Reynolds number between 9.5 x 10^ and 5.2 x 10^. Here the larger Reynolds 
number is close to the value describing the flow near the tip of a helicopter 
main rotor with the same chord and a tip Mach number near M = 0.8. 

Although the operating Reynolds number range of the model airfoil was 
selected to simulate full-scale rotor conditions, transition was not expected 
to occur at the same location as on a production rotor blade. This is because 
the rotor segment in figure 6 had a smooth surface finish since it was machined 
airfoil section. Thus, transition would not occur until the start of the ad- 
verse pressure gradient. In contrast production rotors usually have an abra- 
sion strip covering the leading edge of the profile as depicted by the airfoil 
cross-section shown in figure 7. The discontinuity in the joint between the 
abrasion strip and rotor surface can be sufficient to cause laminar to turbu- 
lent transition of the boundary layer upstream of the minimum static pressure. 

Fortunately, the start of the adverse pressure gradient and the abrasion 
strip joint are located within the first 10 percent of the model airfoil chord. 
This is substantiated by figure 8 which shows the location of the minimum 
static pressure station on the airfoil suction side for the various angles of 
attack tested in the present study. Also shown is the approximate location of 
the abrasion strip joint on the specific production rotor which the model air- 
foil segment simulated. 

To ensure that the transition point always occurred at the same point on 
the model airfoil the boundary layers were tripped at the 10 percent chord 
station on both sides of the airfoil. An additional trip, located at x/c = 

0.08, was needed on the pressure surface side to enhance the surface roughness 
at positive angles of attack. Vortex shedding noise observed during pre- 
liminary tests at positive angles of attack indicated that a single trip would 
not provide a turbulent flow at the trailing edge on the pressure side of the 
airfoil. 
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The boundary layer tripping device employed a thin serrated aluminum tape 
with its serrations facing forward as shown in figure 6. This geometry is 
based on the technique reported by Hama in reference 42. Tape thickness was 
selected to provide minimum surface protrusion to avoid unnaturally large or 
distorted trailing edge turbulent boundary layers. The criteria selected in 
the present study required that the ratio of tape thickness, T, to the local 
boundary layer displacement thickness, 6 , be unity, i.e., T/6 = 1. This 

avoided disturbing the freestream velocity profile outside the boundary layer. 
Local 6 values were estimated using the flat plate boundary layer equation. 
Velocities needed for this calculation were obtained from surface static pres- 
sure curves which were used to calculate the local inviscid velocity field. 

Local Reynolds number based on the velocity at the trip thickness was 
approximately 400 on the airfoil suction side. Hama reported that the serrated 
trip becomes effective at a critical Reynolds number of Re^^ ^ 45. This 
implies that the present trip thickness was sufficient to cause transition. 
Moreover, a constant spanwise transition station was assured because of the 
three-dimensional nature of the boundary layer disturbance (ref. 42). Note 
that the serrated trip approach is more effective than the traditional cylin- 
drical trip device for which 900 (ref. 43) or the distributed roughness 

approach with Re^^ 'v 120 (ref. 43). 

The presence of the above described boundary layer trips on the pressure 
side of the airfoil eliminated the phenomenon of discrete vortex shedding noise. 
Without these trips the boundary layer flow was laminar at positive angles of 
attack. Paterson, et al. (ref. 44) demonstrated that such a laminar boundary 
layer can cause vortex shedding noise at discrete frequencies. 

Since the study of Paterson, et al., vortex shedding has been the subject 
of numerous model rotor, propeller, and airfoil studies. This mechanism, 
however, would not be expected to occur on production main rotors since the 
abrasion strip would cause early transition. Also, the high Reynolds number 
of the flow near the outer portion of the rotor span would result in early 
transition for this dominant noise region of the rotor. Therefore, from a 
practical standpoint vortex shedding is not ah important noise mechanism. 
Consequently, the boundar layer trip on the two-dimensional airfoil ensured 
that the model airfoil tests simulated full-scale operating conditions. Addi- 
tional verification of the absence of vortex shedding was obtained from the 
acoustic spectrum measurements. 


Directional Microphone System 
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The objective of the experimental study was to measure the trailing edge 
noise generated by a two-dimensional segment of a helicopter rotor. The various 
facility background noise sources, shown in figure 9, dominated over the airfoil 



noise. To suppress these extraneous sources and detect only the airfoil noise 
required the use of the directional microphone system illustrated in figure 9. 
This approach was previously shown by Schlinker (ref. 5) to be successful at 
isolating airfoil self noise. A brief discussion will be provided here to 
familiarize the reader with the technique with additional details presented in 
Appendix A. 

Operating Characteristics - The UTRC wide frequency band, highly direc- 
tional microphone system is similar to that previously developed for jet noise 
source location measurements. The experimental technique, based on an optical 
analogy, uses the reflecting surface of a concave mirror to form an image of 
the sound emanating from a segment of the airfoil surface. The present design 
was based on the work by Laufer, Schlinker, and Kaplan in reference 45 (see 
Schlinker (ref. 46) for a more detailed discussion), in which quantitative 
measurements of supersonic jet noise acoustic source distributions were made. 

The 1.067 m aperture spherical reflector, shown schematically in figure 10, 
has a radius of curvature * 1.346 m. A 0.635 cm microphone, situated at the 
focal point, is directed toward the reflecting surface (figs. 11, 12). A 
microphone grid is used in all experiments. The focal point microphone dis- 
tance, r^, from the center of the reflector is controlled by the simple lens 
equation shown in figure 12. 

Test Set-up - The directional microphone test set-up is shown in the 
photographs of figures 5 and 11. The system was traversed parallel to the open 
jet from a station upstream of the airfoil leading edge to a station downstream 
of the trailing edge. Vertical position of the reflector centerline was aimed 
at midspan of the test section. 

To obtain source directivity information in the plane normal to the air- 
foil, the system was aimed at several different angles relative to the open jet 
centerline. As the unit scanned past the airfoil, the output of the focal 
point microphone and the axial position sensor were recorded on magnetic tape. 
Traverse speeds were always less than 0.5 cm/sec. 


Forward Flight Effects 

Shear Layer Refraction and Scattering of Sound - For tests conducted at 
freestream Mach numbers less than 0.1, measurements outside the open jet can 
be used to directly infer the source noise characteristics. However, at 
higher Mach numbers, the open jet technique is influenced by the shear layer 
through which the sound is transmitted. The shear layer refracts, reflects, 
and scatters the sound radiated from the model significantly altering the 
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acoustic source directivity pattern* A brief discussion of the refraction 
corrections is presented here. Changes due to reflection by the shear layer 
were minimal for the angular measurement range investigated here. Additional 
details of the refraction corrections and the resulting ’’apparent” source 
position are given in Appendix B. Scattering effects on directional micro- 
phone measurements are also described. 

An example of the ray path geometry for transmission of sound through a 
shear layer is shown in figure 13. Here, wavefronts propagating in the direc- 
tion, e^, inside the airstream propagate in the direction, 6^, outside the 
open jet. The actual path of the acoustic ray is described by the points SBO. 
The change from 6^ to 0^ is a result of refraction by the shear layer. Without 
flow the sound propagates at angle 0^. 

Directivity information in the present study was obtained over a range of 
angles, 0^. To account for the refraction effect, the directional microphone 
system was rotated in advance to angle, 0^, to obtain the sound initially 
radiated in the direction, 0^. This is illustrated in figure 14(a). Scanning 
the directional microphone past the airfoil trailing edge isolated the sound 
radiated in the direction of interest. 

Airfoil Angle of Attack in Open Jet - Lifting airfoils within an open 
jet induce curvature of the shear layer and deflection of the jet. If the 
wind tunnel has no downstream physical constraint on position of the deflected 
jet, the lifting airfoil rotates the jet until its downward momentum is equal 
to the lift force per unit time. The relative airfoil angle of attack is 
then reduced by the jet deflection. This effect was not accounted for in the 
present study. All angle of attack values, therefore, refer to the geometric 
angle of attack. 


Ins t rumen t a t ion 

Acoustic sound pressure levels at the reflector focal point were measured 
using a 0.635 cm diameter microphone (with grid) at normal incidence. The 
normal orientation introduced a frequency dependent microphone response which 
is typically a function of the angle between the normal to the diaphragm and 
the acoustic ray direction. This effect was included in the directional 
microphone gain calibration which compared the focal point microphone output 
with, an omnidirectional microphone output at grazing incidence. The gain 
measurement, which was previosuly reported in reference 5, was repeated in the 
present study as a check on the directional microphone system (see Appendix B) . 

The focal point microphone signal was amplified and recorded on magnetic 
tape. The frequency response of the FM tape system was flat over the frequency 
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range investigated here. The output from a linear potentiometer, used to 
sense the directional microphone position, was also recorded. The data was 
processed with a 1/3 octave band analyzer. An X-Y plotter displayed the 
selected 1/3 octave band sound pressure level as a function of reflector 
aiming position. 

The slew rate associated with the continuously varying focal point micro- 
phone output was less than the maximum acceptable rate for the 1/3 octave band 
analyzer. This was verified experimentally by stopping the transverse and 
checking for an analyzer output lag. Although this was not a problem in the 
present study, directional microphone measurements of a moving source could 
encounter such difficulties. 

Mean velocity and turbulence intensity measurements at the airfoil trail- 
ing edge were acquired by a 0.025 mm diameter hot film probe operating in 
conjunction with a constant temperature anemometer system. The probe was 
calibrated directly in the tunnel test section. 


Test Program 

Airfoil trailing edge noise measurements were obtained over a range of 
Mach numbers, propagation angles, and airfoil angles of attack to assess the 
dependence on these parameters in the scaling laws and the trailing edge noise 
theory. Table I lists the test conditions. 

The a = -0.4® angle of attack case in Table I corresponds to zero section 
lift on the airfoil. At this operating condition the turbulent boundary 
layers on the pressure and suction side of the airfoil would be expected to 
be approximately the same. In this case, both boundary layers would generate 
the same noise characteristics. 

The a = 12® angle of attack provided the highest section lift without 
stalling the airfoil based on available airfoil lift curves. Mach number was 
limited to M = 0.3 for this operating condition because of the dramatic in- 
crease of facility background noise. This was because lifting airfoils at 
high angles of attack in an open jet induce curvature of the shear layer and 
deflection of the open jet. This deflection impinges the high velocity jet 
on the downstream collector generating sufficient noise to dominate over the 
signal enhancement obtained from the directional microphone. 
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Directivity measurements covered a large angular range for M = 0.3, but 
at higher Mach numbers only limited measurements could be obtained. In the 
upstream quadrant the limitation was due to the large refraction angle change. 
For example, at M = 0.55, sound radiated at 0^ = 90° is refracted to 0^. = 142°. 
Placement of the sideline directional microphone traverse system for measuring 
the sound radiating at 0^ = 142° was not possible due to the presence of the 
upstream anechoic chamber wall. The dependence of 0^ on 0^ and M is given in 
Table II. In the downstream quadrant, measurements at small values of 0^, 
corresponding to small 0^ angles, were limited by the dominance of the open 
jet background noise. Since jet noise varies with velocity to the eighth power 
at high jet velocities, it dominated over the 5th power dependence of the air- 
foil trailing edge noise. 
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AIRFOIL BOUNDARY LAYER CHARACTERISTICS 


Objective 

The objective of the airfoil boundary layer measurements was to document the 
boundary layer thickness, 6 , at the trailing edge of the model airfoil. Direct 
measurements of this parameter were needed as an input to the trailing edge noise 
predictions. In addition, such measurements provided an assessment of the accuracy 
of using flat plate boundary layer calculations as an input to the noise predictions. 
If the flat plate boundary layer characteristics could be shown to predict the mea- 
sured acoustic data with sufficient accuracy, then noise predictions would not re- 
quire detailed velocity profile information at the airfoil trailing edge. This 
could, of course, only be expected to hold at small angles of attack where the flow 
is attached to both surfaces. 

Boundary layer thickness is considered to be the length scale parameter in- 
volved in the trailing edge noise generation process. This follows from previous 
studies which showed that the size of the turbulent structure generating the noise 
is of the order of 6. For example the flow vizualization experiments of Yu and Joshi 
(ref. 9) revealed the presence of large-scale coherent motions in the outer part of 
the turbulent boundary later. The approximate scale of the coherent structure was of 
the order of the boundary layer thickness. Correlation studies concluded that the 
passage of these large-scale eddies over the airfoil trailing edge generated the 
trailing edge noise. 

Previous investigators have used either the boundary layer thickness, 6, or 
the displacement thickness, 6*. The latter parameter is usually chosen when the 
velocity profile shape makes it difficult to identify the edge of the boundary lay- 
er. The displacement thickness however, does not fully represent the physics of 
the noise generation mechanism as indicated by the above discussion. 

Approach 

Although the boundary layer thickness controls the noise generation process, 
it is often a difficult parameter to identify. A comparable length-scale parameter 
which can be easily measured is the location of the turbulent boundary layer inter- 
mit tency interface. The mean interface location is determined by the value of the 
turbulent intermittency parameter, y. This parameter defines the percentage of time 
the signal from a hot film probe displays an intermittent or turbulent signal when 
located near the outer edge of the boundary layer. Since the value of y approaches 
zero at the edge of the boundary layer the location at which the hot film signal 
becomes steady can be used to identify the edge of the turbulent boundary layer. 


Klebanoff (ref. 48) indicated that for a flat plate the height above the 
surface corresponding to an intermittency value of y = 2% occurred at a distance 
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of approximately of 1.07 6. This is illustrated schematically in figure 15 
where the boundary layer thickness, based on U/U^ = 0.99, is compared with the 
intermittency interface location, 6^. Notice that the distance at which the 
hot wire signal is unsteady for 2% of the time is comparable to the boundary 
layer thickness. Since 6^ is easily determined using the visual display of 
the hot film signal on an oscilloscope this method was employed to obtain an 
approximate measurement of the boundary layer thickness in the present study. 

The above argument showing the equivalence of 6 and 6 ^ was previously ver- 
fied only for a flat plate boundary layer. To assess the accuracy of this 
equivalence for an airfoil, mean velocity profiles were also measured in the 
present study. Comparisons of the value of 6 ^ with the normalized velocity 
profile curves indicated that the intermittency interface occurs near U/Uq = 

0.99. Thus, at the airfoil trailing edge in accordance with the flat 

plate data of Klebanoff. The intermittency interface distance, therefore, 
provided a comparable but simplified method for determining the boundary layer 
thickness. 

All velocity profile surveys were conducted at a distance of 0.007 chords 
downstream of the airfoil trailing edge. Admittedly, these measurements corre- 
sponded to wake velocity profile surveys instead of boundary layer profile 
surveys on the airfoil surface. However, the close proximity of the measure- 
ment station to the airfoil trailing edge ensured that the outer edge of the 
wake velocity profile was similar to the boundary layer just upstream of the 
trailing edge. This was expected to be true for the small angle-of-attack 
attached flow conditions investigated in the present study. 

Calibrations of the hot film probe employed in this study were conducted 
in the open jet test section. A commercially available linearizer circuit was 
used to obtain a linear velocity versus voltage response curve over a range 
of 0-191 m/sec. However, for velocities between 137 m/sec and 191 m/sec the 
circuit output became nonlinear. This occurred due to the decrease in static 
pressure inside the open jet test section at high velocities. The pressure 
changes resulted in a decrease of the fluid density which influences the heat 
transfer rate on the hot film. Such changes were not accounted for in the 
analog circuitry of the commercial linearizer unit since most hot film anemometer 
measurements are conducted in constant static pressure environments. For 
example, measurements of the downstream flow field development in a model jet 
are typically conducted at ambient pressure. Fortunately, the effect of 
static pressure changes were included in the present calibration data since 
the hot film probe was calibrated in the open jet test section. Absolute 
velocity profiles on the airfoil surface were then determined from the measured 
calibration curve. It is recognized that due to the nonlinear response of the 
anemometry circuitry at high Mach numbers turbulence measurements were only 
meaningful for m50.4. 
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Results 


Mean Velocity Profiles and the Intermittency Interface Location - Figure 16 
shows the measured mean velocity profile and the location of the intermittency 
interface on both sides of the model airfoil at a = -0.4°. All dimensions have 
been normalized by the airfoil chord and velocities have been normalized by the 
maximum measured freestream velocity. Velocity profile data near y/c = 0 was 
not needed for the assessment of the noise predictions in the present study. 
However, such information was essential to determining the position of the hot 
film probe relative to the airfoil surface. Since velocity profile measurements 
were made at a minute distance downstream of the trailing edge, it was necessary 
to reference all transverse distances to a hypothetical extension of the airfoil 
chord. Thus, the y/c = 0 station in figure 16 coincided with the hypothetical 
airfoil extension. 

The lack of a physical reference point downstream of the trailing edge 
made it difficult to identify the location of airfoil chord extension. This 
could only be determined visually be requiring the hot film to bisect the pro- 
jected airfoil surface area when viewing the airfoil trailing edge from the 
downstream end. To provide a check on this approach the measured velocity pro- 
files near y/c = 0 were used to identify the trailing edge stagnation point. 

Two straight lines were drawn through the suction and pressure side data points 
near y/c = 0 in figure 16. The intersection of these lines approximated the 
location of the stagnation point at the airfoil trailing edge. For the zero 
angle-of -attack case shown in figure 16 the location of this stagnation point 
agreed closely with the y/c = 0 reference point. Measurements of 6 ^ were, 
therefore, obtained relative to the desired extension of the airfoil chord. 

Comparing the location of the intermittency interface with the velocity 
profile curves near U/U^ = 0.99 indicates that 6 ^ provides a comparable measure- 
ment of the boundary layer thickness. The replacement of 6 by 6 ^ has, there- 
fore, been verified for the a = -0.4° angle of attack test condition. 

Also shown in figure 16 are the normalized turbulence intensity profiles 
defined by As expected from the criterion for selecting, the 

fluctuating measurements approached zero at the intermittency interface. This 
once again confirmed that identifies the edge of the turbulent boundary layer. 
It should be noted that turbulent fluctuation measurements at the higher Mach 
numbers of M = 0.43 a nd 0.5 c orre sponded to the turbulent mass fluctuations 
defined by the ratio (P )/(P )inax* Although this ratio is plotted in 
figures 16(c) and 16(d) absolute values could not be specified due to the non- 
linear characteristics of the hot film response at the high velocities. 

A comparison of 6 ^ on both airfoil surfaces in figure 16 indicates that the 
boundary layers are similar at a = -0.4°. This is demonstrated in figure 17(a) 
where the values of S^/c are plotted for the pressure and suction side surfaces 
at the various Mach numbers investigated in the present study. Also shown in 
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figure 17(a) is the flat plate boundary layer thickness calculation for 6/c. 

Here the turbulent boundary layer is assumed to start at the transition point 
determined by the boundary layer trip located at x/c = 0.1. The good agree- 
ment between -the predicted values of 6/c and the measured values of <5^/c 
indicates the applicability of the flat plate prediction at zero lift condi- 
tion when both boundary layers develop under the influence of approximately 
the same adverse pressure gradient. 

The effect of Mach number or Reynolds number is weak as shown by the 
experimental data in figure 17(a). This agrees with the flat plate turbulent 
boundary layer prediction in which 6 varies according to 

Figure 18 shows the boundary layer velocity profiles for a = 7.6° angle 
of attack. The boundary layers are no longer symmetric as in the a = -0.4° 
case. Here the airfoil trailing edge stagnation point, approximated by the 
intersection of the two straight lines, shows a small disagreement with the 
y/c = 0 station determined during the experiment. The offset was used to 
correct the measured values of 6^ to account for the apparent error in deter- 
mining the y/c = 0 reference station. The advantage of measuring the complete 
velocity profile is now evident. 

The mean velocity profiles in figure l8 indicate that the boundary layers 
are no longer S57mmetric. The intermittency interface station, however, still 
occurs where the mean velocity profile as 3 anptotically approaches U/U^ = 0.99. 
Thus, 6^, continues to provide a representative measurement of the boundary 
layer thickness, 6. 

The variation of 6^^ with Mach number is summarized in figure 17(b) for 
the a = 7.6 case. Similar to the a = -0.4°, changes with increasing Mach number 
are small. The flat plate boundary layer prediction would, obviously, be unable 
to predict the upper and lower surface values of 6^. 

After demonstrating that 6 can be replaced by 6^ at a = -0.4° and 7.6° 
measurements at a = 12° were limited to determining the intermittency inter- 
face location. Figure 17(c) shows the experimental results at this highest 
angle of attack. As expected the suction surface boundary layer thickness 
has increased significantly relative to the a = -0.4° case while the pressure 
side boundary layer thickness has decreased. 

Acoustic Scaling Law Dependence on Boundary Layer Thickness - Existing 
predictions for airfoil self noise contain a single boundary layer thickness 
parameter for each side of the airfoil. Unfortunately, at high-lift conditions 
the boundary layers on both sides of the airfoil are significantly different. 

One approach to circumventing the resulting dilemma is to use the average of 
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the suction and pressure side boundary layer thickness values. The Mach 
number and angle-of -attack dependence of this average value was obtained 
for the 6^ values in figure 17 and plotted in figure 19. A mean line has 
been drawn through the measured data points for use in later data reduction. 

For a = -0.4® and 7.6° the average value of 6^ is within 20% of the calculated 
flat plate boundary layer thickness when M f 0.3. Note that the effect of angle 
of attack on the radiated noise is considered to be weak for small a. This is 
because the boundary layer thickness parameter appears as a logarithmic depen- 
dence in the trailing edge noise scaling laws so that 20% error in 6^^ introduces 
at most a 0.8 dB error in the amplitude scaling. 

At the highest angle of attack in figure 19 the average value of 6^ no longer 
agrees with the flat plate boundary layer calculation for 6. Differences between 
data at the two lower angles of attack and this highest angle of attack are 
presently not understood. 


35 



DATA ACQUISITION AND DATA REDUCTION PROCEDURE 


A demonstrated experimental approach to isolating turbulent boundary layer 
noise is through use of two-dimensional airfoil sections tested in an acoustic 
wind tunnel with low turbulence levels. In the present study a directional 
microphone was needed to separate the turbulent boundary layer noise from other 
facility background noise sources. Directivity information was obtained by 
orienting the directional microphone system at different viewing angles relative 
to the airfoil. Examples of the data acquired from the directional microphone 
system are given below. Also, a brief discussion of the data reduction procedure 
is provided to describe how the source distribution measurements are converted 
to absolute sound pressure levels associated with trailing edge noise. 


Directional Microphone Data Acquisition Procedure 

For each Mach number and directivity angle 0^, there exists a location in 
the open jet corresponding to the apparent source position. This apparent 
source location controls the sideline traverse distance in figure 13 since 
the reflector-to-apparent source distance was required to be R = 2.07 m or 2.81 m. 
Once the sideline distance was determined the reflector centerline orientation 
relative to the open jet axis was adjusted to coincide with the transmission 
angle, 0^, determined from refraction calculations. As shown in Table II, 0^ 
varies with 0^ and M. 

After selecting the reflector sideline distance and orientation angle the 
directional microphone system was traversed past the isolated airfoil. Figure 20 
is an example of the measured source distribution as the directional microphone 
system was traversed parallel to the open jet centerline. Such curves were ob- 
tained at each 1/3 octave band frequency between 1 kHz and 50 kHz. Each trace 
was determined by playing back the analog record from the tape recorder and 
passing it through the selected 1/3 octave band. 

The first peak near the tunnel exit plane in each trace represents the open 
jet shear layer noise and nozzle lip noise. The spherical directional microphone 
system is sensitive to these acoustic sources because of the finite depth of 
field. Note that a large distance separated the nozzle exit from the airfoil 
trailing edge to avoid contamination of the airfoil source distribution by the 
tunnel lip noise. 

The second peak located downstream of the airfoil trailing edge location in 
figure 20 is associated with acoustic radiation from the airfoil model. The dis- 
placement between the trailing edge and the apparent source position agrees with 
the expected shift due to convection of sound waves by the tunnel freestream. 

The downstream displacement of the apparent acoustic source position relative 


36 



to the airfoil trailing edge was determined analytically using the simple 
expression Ax = h^U/c^^. Here h^^is the trailing edge-to-shear layer 
distance (see fig. 44) for a ray propagating at angle 0^ when M = 0. The 
parameter, Ax, represents the distance which the reflector, initially aimed at 
the trailing edge, must be moved downstream to measure the noise radiation in 
the presence of flow (fig. 44). Comparisons showed good agreement between the 
measured and predicted values of Ax. The good agreement between theory and 
experiment also provided an indirect confirmation of the refraction prediction 
for the angle Sf 

To verify that the airfoil trailing edge noise mechanism was responsible 
for the second peak in each distribution, the airfoil was removed and the tra- 
verse was repeated for a few selected cases. The f - 5 kHz curve in figure 20(a) 
shows the resulting source distribution which is referred to as the facility back- 
ground noise. Lip noise, open jet shear layer noise, and collector noise are 
the only contributors to this background noise. 

One advantage of the directional microphone is now evident. The reflector, 
with its large gain, enhances sound arriving from a source on the mirror axis 
while rejecting sound arriving from other off-axis sources. In comparison, a 
single omnidirectional microphone would measure all three background noise 
sources in addition to the trailing edge noise. As previous investigators 
have indicated, the trailing edge noise source would be masked by the facility 
noise in this case. 


Calculating Trailing Edge Noise Spectra 
from the Directional Microphone Data 

Acoustic source distributions like those in figure 20 are typical of the 
features observed over the range of test conditions. From such curves the far- 
field trailing edge noise spectra were calculated. Specific equations developed 
for the data reduction purposes were presented by Schlinker in reference 5. 

The calculation procedure uses the absolute sound pressure level measured at the 
peak of the source distribution associated with the trailing edge noise radiation. 
Due to the presence of background noise sources it is necessary to correct the 
measured sound pressure level. Ideally, this would be obtained by subtracting the 
value of the background noise curve directly under the peak associated with the 
trailing edge noise. 

Unfortunately, the background noise curve shown in figure 20(a) is a function 
of the airfoil angle of attack. This is because for high values of lift coeffi- 
cient the airfoil deflects the open jet sufficiently to cause increased collector 
impingement noise. Removal of the airfoil to determine the open jet background 
noise levels is, therefore, possible only at 2ero lift conditions where the air- 
foil has no affect on the impingement noise. For this reason, a direct measure- 
ment of the facility background noise curve associated with finite lift conditions 
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was not possible. Instead, straight lines were drawn connecting the broadband 
noise levels on either side of the trailing edge noise peak in figure 20(b). The 
assumed sound pressure level beneath the peak represented the facility background 
noise. Based on the flat shape of the curve on either side of the peak in 
figure 20(b) , a straight line was adequate for defining the facility background 
noise distribution. 

The directional microphone data reduction procedure described in reference 5. 
determines the equivalent trailing edge noise spectrum which would be measured by 
a hypothetical omnidirectional microphone situated at the reflector station at 
which the peak occurred in figure 20. The equivalent noise spectrum in this case 
represents the sound field sensed by an omnidirectional microphone situated at 
the observer position (0) in figure 13. Applying the refraction angle and 
amplitude corrections (Table II) to the resulting experimental data converts the 
measurements at station 0 to the sound pressure level that an omnidirectional 
sideline microphone would sense at station D in figure 13. These corrections 
transform the experimental data to a measurement system in which the microphone 
is located inside the flow field and fixed relative to the noise source. Thus, 
the converted data corresponds to an airframe flyover measurement in which the 
microphone, situated at angle 6^, is moving at the aircraft flight speed in a 
direction parallel to the flight path. This coordinate system is often referred 
to as the present coordinate system. 

To facilitate comparisons with theoretical directivity predictions, it is 
convenient to convert all present coordinate system measurements to a retarded- 
time coordinate system. This latter coordinate system effectively moves the 
noise source downstream to its relative physical location in addition to changing 
the radiation angles from 6^ to 6^. The retarded source coordinate defines the 
position of the source when it emitted the sound sensed at the microphone station. 
Thus, conversion of the data to this coordinate system would represent the 
sound pressure level measured with the microphone fixed and the source moving 
at speed, M, in figure 13. The Galilean transformation which results in the 
source moving also results in a zero open jet velocity. 

The ray path angles linking the present and retarded coordinates are given 
by the frequency independent equation 

cos = cos dc sin^ - m sin^ 6^ (60) 

In addition to converting to the retarded source angles a sound pressure level 
correction was applied to the data. Measured sound pressure levels were converted 
from the sideline location (point D in fig. 13) to a common source radius in the 
retarded coordinate system given by r^ * 3 m. The correction was determined from 
the equation 
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( 61 ) 


A Doppler shift was not included in the transformed data since the theoreti- 
cal predictions of trailing edge noise for the isolated airfoil test conditions 
did not include this effect. The resulting angle and amplitude changes between 
the measured data and the retarded frame of reference are summarized in Table II. 

After correcting for refraction effects, converting to retarded source 
coordinates, and correcting for the directional microphone resolution and gain 
characteristics, the 1/3 octave band sound pressure levels associated with each 
test condition were plotted as a function of Strouhal number. In this case the 
Strouhal number was based on the parameter St = f 6 /Uq. Here, the length scale 
parameter, 6, is determined from 6 = 6j^/1.07 where 6^ corresponds to the average 
of the measured intermittency interface distance presented previously in 
figure 19. The frequency, f, corresponds to the 1/3 octave band center fre- 
quency and the velocity, U, corresponds to the freestream velocity of the open 
jet test section. 
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EXPERIMENTAL ASSESSMENT OF TRAILING EDGE NOISE RADIATION 
FROM A TWO-DIMENSIONAL AIRFOIL SECTION 


Objective 

The objective of this phase of the study was to experimentally assess the 
noise generated by a turbulent boundary layer flowing over the trailing edge 
of an airfoil. Tests conducted over a range of Mach numbers, airfoil angles 
of attack, and directivity angles were used to assess the dependence of the 
trailing edge noise on these physical parameters. Conclusions obtained from 
this assessment were used to develop a scaling law based on the different 
physical parameters. Included in the data base used to assess the scaling 
law were the trailing edge noise measurements reported by previous investigators 
in references 5 and 10. These different experimental studies provided data 
with widely varying freestream Mach number, airfoil chord, and airfoil pro- 
file values, thereby generalizing the scaling law applicability. The experi- 
mental measurements cited in references 5 and 10 were also compared to the 
fundamental analytical approach described in the section titled Theoretical 
Formulation of the Trailing Edge Noise Mechanism. This provided a preliminary 
evaluation of the accuracy of the theoretical prediction procedure. A detailed 
assessment is presented in the following section titled Trailing Edge Noise 
Predictions. 


Angle-of-Attack Assessment 

Approach - The effect of angle of attack on airfoil trailing edge noise 
was assessed by comparing acoustic spectra for different values of a. Compari- 
sons were obtained over a range of Mach numbers as angle of attack varied from 
a - -0.4® to 12®. All spectra corresponded to a fixed radiation angle of 
0^ = 70®. When transformed to retarded coordinates, the radiation angle, 0^, 
varied with Mach number according to equation (60). Angle 0^ is indicated in 
each data set. Note that the spectra have been converted to the acoustic 
radiation from a unit length (1 ft or 0.33 m) of the airfoil. 

Prior to evaluating the experimental results, a brief discussion of the 
spectrum measurement range is needed. Due to the lack of data points at low 
Strouhal numbers, it was not possible to identify the trailing edge noise spec- 
trum peak (see, for example, fig. 21). This limitation occurred because of 
the low signal-to-noise ratio for directional microphone measurements at fre- 
quencies below 1.25 kHz. In this range, the reflector gain was too small to 
raise the trailing edge noise level above the facility background noise which 
radiates directly to the focal point microphone. Also, the large depth of 
field at low frequencies resulted in the reflector sensing the open jet shear 
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layer noise. Despite these limitations to measuring the complete acoustic 
spectrum, the experimental data obtained from the present study was sufficient 
to isolate the parametric dependence of trailing edge noise. This will become 
apparent in the following discussions. Also plotting the data in a normalized 
format, as is done in the section titled Generalized Scaling Law, shows that 
the lowest measured Strouhal numbers approached the spectrum peak in most 
cases. 

Experimental Results - The general trend in figure 21 shows an increasing 
sound pressure level in the low Strouhal number range as the airfoil angle of 
attack increases. This dependence can be explained by the change in the 
boundary layer thickness on the suction side of the airfoil. As angle of 
attack increases, the adverse pressure gradient on the suction side creates a 
thicker boundary layer (see, for example, fig. 17). Since turbulent eddy size 
scales with the boundary layer thickness, increasing a creates larger eddies 
on the suction side of the airfoil. The presence of these large eddies results 
in a boundary layer pressure spectrum dominated by low frequencies which in 
turn increased the low frequency acoustic radiation. 

The above physical explanation is substantiated by the flat plate turbu- 
lent boundary layer wall pressure fluctuation study by Schloemer (ref. 49). 

The tests were conducted with zero, adverse, and favorable pressure gradients. 

The results indicated that the effect of an adverse pressure gradient on the 
nondimensionalized wall pressure spectral density is to increase the low fre- 
quency content without influencing the high frequency portion appreciably. In 
contrast, a favorable pressure gradient creates a sharp decrease in the spec- 
trum at high frequencies. Since the turbulent boundary layer pressure spectrum 
controls the acoustic radiation, similar trends would be expected in the acous- 
tic far field. 

Indeed, the data in figure 21 demonstrates an Increase in low frequency 
acoustic energy, thereby, verifying the postulated link between boundary layer 
thickness and spectrum content. At high frequencies the expected decrease in 
sound pressure level did not occur. Instead, in the mid-frequency range, the 
noise levels remained constant. One possible explanation is that changes in 
the suction side and pressure side boundary layers compensated each other in 
the mid-frequency range, given by St ^ 0.8. For even higher Strouhal numbers, 
a different explanation is, however, needed. Now an increase in angle of attack 
also increases the noise radiation (see, for example, figs. 21(b) and 21(c)). 

These changes, which occurred mainly for M - 0.3, are presently not understood. 
Fortunately, the high frequency portion of the spectrum generated at the full- 
scale Reynolds number operating conditions tested here, would not be audible 
(f ^ 20 kHz). Thus, it was not critical to understand the observed changes. 

These changes were, however, considered to be due to a different source mechan- 
ism since the spectrum changes were centered near St = 3 as shown by figure 29(c). 
This Strouhal value is approximately a factor of 10 above the value associated 
with the trailing edge noise mechanism. 
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One noise mechanism which was evaluated as a possible explanation for the 
high frequency anomaly in figure 29(c) was vortex shedding due to the blunt- 
ness of the airfoil trailing edge. Assigning a Strouhal number of 0.2 to this 
mechanism and using the source frequency associated with the center frequency 
of the hump (f = 31.5 kHz in fig. 21(c)) provided a calculated dimension of 
0.065 cm. This dimension is of the same order of magnitude as the airfoil 
trailing edge thickness which is 0.038 cm. Blunt body vortex shedding is 
therefore, a possible mechanism responsible for the high frequency portion of 
the spectra in figure 21. 


Directivity Assessment 

Approach - Directivity patterns were evaluated in the overhead plane at 
specific Strouhal numbers for various angles of attack and open jet Mach num- 
bers. Sound pressure levels were normalized by the level measured at a speci- 
fied retarded angle. The value of 0^ selected to normalize the data was close 
to 90° in each case. Theoretical directivity curves were also normalized by 
the same reference angle. 

Evaluating the directivity pattern for various Strouhal numbers provided 
several advantages compared to the standard method of using the overall sound 
pressure level (OASPL) . First, it permitted evaluating the directivity depen- 
dence on Strouhal number. Second, it permitted selecting the segment of the 
spectrum responsible for trailing edge noise while rejecting the noise postulated 
to be due to the vortex shedding source mechanism described in the previous 
subsection. Finally, the directivity dependence could be assessed despite the 
inability to determine OASPL values due to the limitations in the spectrum 
measurement range. 

Directivity measurements were referenced to the open jet centerline which 
was within -0.4° of the airfoil chord line at zero lift conditions. The small 
changes in airfoil orientation due to Increasing angle of attack were not in- 
cluded in the data reduction. The resulting error was small due to the slow 
variation in the radiation field with 0^ for the angular range Investigated 
here. At most, the error was equivalent to 0.75 dB at the highest Mach number 
of M * 0.5. 


Experimental Results and Comparison with Present Theory - Data showing the 
directivity pattern for various Strouhal numbers are shown in figure 22 for a 
range of velocities with a = -0.4°. The measurements show little difference in 
the directivity shape for Strouhal numbers below St = 1.0. The similarity of 
the data curves indicates that the directivity pattern is independent of 
Strouhal number. 
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The effect of angle of attack is determined by comparing the a = 0.4® case 
in figure 22(b) with the a * 7.6® case in figure 23. For the limited angular 
range available in figure 23, the directivity pattern is approximately similar 
to the measurement at zero lift coefficient. This constant shape of the radia- 
tion field indicates that the directivity pattern is independent of angle of 
attack. 

Also shown in figures 22 and 23 is the analytical trailing edge noise 
directivity pattern determined from equation (32) of the present study. This 
radiation pattern, which is independent of Strouhal number, was normalized to 
the same retarded angle as used in the experimental data analysis. It is seen 
that the measured sound field agrees with the predicted shape to within 1 dB 
in most cases. 


Mach Number Assessment 

Approach - Mach number scaling laws were assessed experimentally using 1/3 
octave band spectra like those shown in figure 24. Sound pressure levels 
measured at a fixed Strouhal number were plotted as a function of increasing 
Mach number. The influence of angle of attack was also evaluated by comparing 
the Mach number dependence at a ® -0.4° and 7.6°. At each Mach number the 
sound pressure level plotted for a specified Strouhal number represented the 
1/3 octave band measurement at the Strouhal number closest to the specified 
value. This occurred because the fixed 1/3 octave bands in the spectrum 
analyzer could not be expected to be centered at the precise frequency calcu- 
lated from the Mach number and boundary layer thickness using the equation 
f = SfU/6^. 

Since the directivity dependence was evaluated in the previous section, 
the present discussion will be limited to the Mach number dependence at a sin- 
gle retarded angle of 6^ = 98®. This was considered close to 6g = 90° where 
the directivity function is independent -of Mach number. The value of 6^ * 98° 
was selected based on the availability of data at this angle for the Mach num- 
ber of Mg “ 0.5. Because measured data at lower Mach numbers was obtained at 
other angles (see Table II for 6^ values), it is necessary to introduce direc- 
tivity corrections to convert the lower Mach number data to 0g = 98°. All 
corrections were based on the theoretical directivity patterns. 

Experimental Results - Figure 24(a) shows the variation of sound pressure 
level versus Mach number for a range of frequencies with angle of attack set 
at a “ -0.4°. Within the accuracy of the figure, the general trend shows good 
agreement with a velocity scaling law of M^ (see solid lines in fig. 24(a)) 
for Strouhal numbers varying over a factor of approximately 10. Note that low 
Strouhal number data for the case of M^ * 0.1 was not available due to the fre- 
quency limitations of the original data. 
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Increasing the airfoil angle of attack to a = 7. 6® in figure 24(b) 
continued to give good agreement with the scaling law. From this it can 
be concluded that at each angle of attack the Mach number dependence is the 
same. Of course, the absolute value of the measured sound pressure level at 
a fixed Strouhal number continues to be a function of both a and M. This 
dependence is demonstrated by the change in spectrum shape and absolute 
amplitude as angle of attack and Mach number varied in figure 21. 

The above described M^ sound pressure level dependence agrees with the 
theoretical analysis of the present study. This is verified in the section 
titled Trailing Edge Noise Predictions. 


Comparison with Previous Investigators 

Experimental and theoretical results of the present study can be compared 
with the results of previous investigators. The following subsection concen- 
trates on the angle of attack, directivity, and Mach number dependence. Com- 
parison of the measured acoustic spectra with spectra reported by previous 
investigators are presented in the subsection titled Spectrum Shape Assessment. 
Finally, comparisons with existing scaling laws are given in the subsection 
titled Generalized Scaling Law. 

Angle of Attack Dependence - The observed increase in sound pressure level 
at low Strouhal numbers as angle of attack increases agrees with the trends 
reported earlier by Schlinker in reference 5. The study by Brooks and Hodgson 
(ref. 10) also showed that angle of attack changes affect primarily the low 
Strouhal portion of the spectrum. Their results, however, showed a tendency 
for the spectrxim peak to shift to lower frequencies in addition to experiencing 
a small decrease in the 1/3 octave band sound pressure level. 

Brooks and Hodgson also observed a distinct peak in their acoustic spectra 
at high frequencies similar to the trend in figure 21(c). Their experiment 
was able to demonstrate that this additive contribution to the spectrum was 
generated by vortex shedding due to the airfoil trailing edge bluntness. By 
increasing the ratio of the trailing edge thickness, T, to boundary layer dis- 
placement thickness, 6*, from T/6* = 0 and T/6* = 0.28, the hump was seen to 
increase in amplitude. 

Directivity Dependence - The analytical directivity pattern described by 
equation (32) of the present study can be compared with the shape predicted 
by the theory of Howe (ref. 11). If the Kutta condition is applied at the 
airfoil trailing edge and the wake convection velocity parameter is assumed 
equal to the turbulent boundary layer convection velocity, the theory of Howe 
gives the same directivity shape as equation (32). If, however, the wake 
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convection velocity in the theory of Howe is considered to be zero, as used by 
Brooks and Hodgson (eq. (31), ref. 10), then the analytical equation describing 
the directivity pattern is different from radiation field defined by equation 
(32) of the present study. 

Although not specially formulated as theoretical study of airfoil trailing 
edge noise, the recent work of Goldstein (ref. 33) also provides directivity 
predictions. The objective of the study was to predict the radiation pattern 
associated with a turbulent shear flow at the edge of a semi- inf inite plate. 

To obtain directivity information relevant to the present study required re- 
moval of the refraction affects which were included in the analysis of ref- 
erence 33. Correcting for shear layer refraction using the approach developed 
by Amiet (ref. 34) resulted in an analytical expression identical to equation 
(32) of the present study. The similarity with the analysis of Goldstein 
occurred only for the high frequency (or high Strouhal number) solution. At 
low frequencies (low Strouhal number) , wehre the acoustic wavelength in the 
analysis of Goldstein was equivalent to the transverse dimension of the shear 
flow, internal reflections by the shear layer complicated efforts to correct 
for the presence of the shear layer. In this case, the analysis of Goldstein 
could not be modified to permit comparison with the analysis in the present 
study. 

With several different theoretical directivity patterns available in the 
literature, it is worthwhile to evaluate further the differences between them. 

A comparison of the directivity pattern obtained from equation (32) and the 
prediction of Brooks and Hodgson (eq. (31), ref. 10) is given in figure 25. 

The parameter, 5(0^, r/2), represents the directivity pattern variation with 
angle 0^. In both cases the calculation for D(0g, ii/2) assumed a turbulent 
boundary layer convection velocity to freestream velocity ratio of 0.6. At a 
Mach number of = 0.2, corresponding to the maximum operating condition 
tested in reference 10, figure 25(a), shows that the two analyses (curves 1 and 
2) give almost identical results. Yet, the two analyses use different values 
for the wake convection velocity. The present study assumes that the wake 
velocity is equivalent to the freestream convection velocity while Brooks and 
Hodgson assign a value of zero to the wake velocity. Thus, good agreement 
between the predicted and experimental directivity patterns is insufficient 
to assess the value selected for the wake convection velocity. 

At a higher Mach number of Mq = 0.43, corresponding to one of the test 
conditions evaluated in the present study, the source directivity patterns 
described by curves 1 and 2 show a different shape (fig. 25(b)). When the Mach 
number is increased further to M^ = 0.86, the differences between curves 1 and 
2 becomes even larger, especially in the forward arc (fig. 25(c)). For example, 
at 0g = 140®, a difference of 10 dB exists. The above described differences 
between the present analysis and the prediction of Brooks and Hodgson 
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demonstrates the need for experimental assessment of the trailing edge noise 
directivity pattern at high Mach numbers. Direct verification using isolated 
airfoils is presently not available. 

Mach Number Dependence - Mach number to the fifth power scaling laws were 
obtained by several previous investigators. Brooks and Hodgson (ref. 10) 
measured an OASPL dependence of Schlinker (ref. 5) also observed a 

fifth power dependence for the maximum sound pressure level in the acoustic 
spectrum. Both of these studies employed isolated airfoils. A similar velocity 
dependence was observed in previous investigations of airframe noise. Fink 
(refs. 6 and 7) reported a fifth power velocity dependence for aerodynamically 
clean conventional low speed aircraft, jet aircraft, and sailplanes. 


Spectrum Shape Assessment 

Objectives - The first objective of this phase of the study was to 
determine if there exists a universal spectrum shape defining the trailing 
edge noise radiation. A second objective was to identify the Strouhal number 
associated with the spectrum peak. Recall that the present study was unable 
to measure the low frequency portion of the spectrum so that the spectrum peak 
could not be isolated. This limitation was overcome by including in the spec- 
trum shape evaluation data sources which clearly identified the spectrum peak. 
These data included the previous reported isolated airfoil measurements of 
Schlinker (ref. 5) and Brooks and Hodgson (ref. 10). If it could be demon- 
strated that these different measurements collapse onto a single curve, then 
the Strouhal number of the spectrum peak could be determined. 

Approach - Acoustic spectra presented here are limited to a single retarded 
angle of 0^ near 90®. In addition, spectrum measurements will be limited to 
a 6® angle-of-attack results. This is because helicopter rotors operate at 
a finite angle of attack. Zero lift coefficient measurements were discussed 
earlier only to provide a complete evaluation of the trailing edge noise 
dependence on angle-of-attack. 

Frequency information in this section of the report continues to be 
presented in terms of a nondimens ional Strouhal number. As discussed in the 
section titled Airfoil Boundary Layer Characteristics, 6 is the appropriate 
length scale to use in the Strouhal number parameter. Using measurements 
reported by previous investigators to check for a universal spectrum shape, 
therefore, required estimating 6 for each study. This was necessary because 
the data reported by Schlinker (ref. 5) did not provide direct measurements of 
the boundary layer thickness. Brooks and Hodgson (ref. 10) presented their 
data in terms of a Strouhal number based on the displacement thickness which 
the authors determined from measured velocity profiles. A direct measurement 
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of 6 was not available possibly due to the asymptotic shape of their velocity 
profile near the edge of the boundary layer. Detailed discussions of the 
procedures used to predict 6 in each of the above studies are given in 
Appendix D, Tables III and IV list the average values of 6 calculated in each 
case. 


Figure 26 shows the measured trailing edge noise sound pressure levels 
obtained from the present study for a = 7.6® angle of attack. The spectra 
were measured at different retarded angles but were not converted to a common 
retarded angle. This was, however, done when the spectra were normalized to 
check for the existence of a scaling law. Note that all acoustic spectra 
represent the acoustic radiation from a unit length (1 ft or 0.3 m) of the 
airfoil. 

Figures 27 and 28 present the acoustic spectra measured by previous 
investigators. Frequency information in each case is given in terms of a 
Strouhal number based on the boundary layer thickness values listed in Tables 
III and IV. These spectra represent trailing edge noise measurements from air- 
foil geometries and chord lengths different from the present study. The follow- 
ing discussion develops the scaling law and tests it against the different data 
sources. 

Scaling Law Development - Dependence of trailing edge noise on Mach number 
and directivity were assessed in the previous subsections. Comparisons between 
theory and experiment indicated that the velocity dependence can be modeled as 
or as U for incompressible flow. In addition, the noise directivity 
pattern can be represented by equation (32) in the plane normal to the airfoil 
model. However, the influence of airfoil profile shape and angle of attack 
have not yet been quantified. These parameters are considered to influence the 
noise indirectly by altering the turbulent boundary layer thickness at the air- 
foil trailing edge. 

Noise dependence on boundary layer thickness was modeled as a linear func- 
tion of 6. A linear dependence on 6 also follows from the theoretical trail- 
ing edge noise formulation. Since 6 is proportional to 6 for the flat plate 
geometry used in the theoretical derivation, 6 can be replaced by 6. 

The theoretical dependence on 6 will become more explicit in the next sec- 
tion titled Trailing Edge Noise Predictions. There the high frequency limit 
of equation (21) will be derived to show the dependence on 6 as well as the 
other physical parameters. The resulting expression guided the selection of 
the specific parametric dependence for use in the scaling law. It should be 
noted that the linear dependence on 6 agrees with the derivation of 
Ffowcs Williams and Hall (ref. 16) in addition to the analysis of Howe (eq. (7), 
ref. 11). 
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Based on the above arguments for the parametric dependence on M, D and 6, 
the overall sound pressure level for a two-dimensional isolated airfoil section 
is given by 


p2 0Cu®15 

*e 


(62) 


Here the parameter s defines the airfoil span and r^ represents the acoustic 
source-to-microphone distance in the retarded coordinate system. In logarith- 
mic form the above equation becomes 


OASPL = 50 log,o 1^+10 log,o % + lO log,o 5 + K, , 

U in knots 


(63) 


where is the unknown constant of proportionality. A reference velocity of 
100 knots was selected for the velocity scaling term to permit comparing the 
value of obtained from the present study with the value previously deter- 
mined in the airframe noise study by Fink (refs. 6andT). For this reason the 
velocity in equation (63) must be in terms of knots. Note that equations (62) 
and (63) apply to an airfoil with zero sweepback angle at the trailing edge. 

If a universal spectrum exists then the OASPL and 1/3 octave band sound 
pressure level are related by a Strouhal dependent function where 


SPL,/3 = OASPL + F, (St) 


(64) 


Here, the parameter SPL represents the 1/3 octave band sound pressure level in 
dB and Fj^ describes the spectrum shape. Solving equations (63) and (64) for 
the SPL dependence gives 

SPL^, = 50log,(, + Ologio ^ + lO iog,o B + F,(st) + K, , 

* (65) 

u in knots 

Recognizing that D = 1 at 6g •» 90®, and isolating both Fj^(St) and gives 


F,(St) + K, = SPl ,/3 -50 log,o TO - 10 log,o 
U in knots 


( 66 ) 
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By subtracting the second and third terms on the right side of equation 
(66) from the different 1/3 octave band spectra in figures 26 to 28, the exis- 
tence of a universal spectrum shape can be tested. If the resulting curves 
of Fj^(St) + for each data set collapse on a single curve then the scaling 
laws in equations (64) and (65) will have been verified. Furthermore, the 
existence of a normalized spectrum shape will have been confirmed. 

Normalized Spectra - Figure 29 shows the normalized acoustic spectra ob- 
tained from the present study. The normalization includes conversion of the 
spectra measured at different retarded angles in figure 26 to a common angle 
of 0^ = 90® where the directivity function, 5 , is unity. Theoretical direc- 
tivity calculations were used for the data conversion. 

Trailing edge noise spectra in figure 29 tend to collapse onto a single 
curve when normalized by the parameters on the right side of equation (66). 
Although a definite spectrum peak cannot be determined from the data shown 
here, the curve suggests that the peak is near St = 0.1. This is similar to 
the value of St = 0.1 determined by Fink (refs. 6 and 7) in the study of 
trailing edge noise generated by fullscale clean airframe configurations. 

Normalized acoustic spectra obtained from the earlier data reported by 
Schl inker are shown in figure 30. The data collapses onto a single curve with 
the spectrum peak occurring near St = 0.1. Values of the parameters used in 
the normalization are given in Table III. 

Figure 31 shows the normalized acoustic spectra obtained from the data 
reported by Brooks and Hodgson for a = 5°. The spectrum peak in figure 31 
occurs near St = 0.14 which is close to the Strouhal peak suggested by 
figures 29 and 30. 

A table of the peak Strouhal numbers associated with the different data 
sources described above gives 


Investigators 


Peak St Value Based on 6 


Present Study 0.1 

Brooks and Hodgson 0.14 

Schlinker 0.1 


Based on the tabulated values the Strouhal peak associated with trailing edge 
noise lies between St = 0.1 and 0.14. This is consistent with the value of 
St = 0.1 obtained from the clean airframe trailing edge noise studies of Fink 
(refs. 6 and 7) but differs from the value reported by Brooks and Hodgson 


49 


(ref. 10). Their study based the Strouhal number on the measured value of 
6* instead of 6 as in the above tabulation. Their range of geak Strouhal num- 
ber values was close to the values tabulated above with St (6 ) varying from 

* ic 

0.07 to 0.1. Here St(6 ) designates a Strouhal value based on 6 . It is not 
clear why the Strouhal values would be similar when 6 is significantly smaller. 
It is possible that the value of 6 was too large as discussed in Appendix D. 

The measurements of Brooks and Hodgson can also be compared with the peak 
Strouhal value (based on 6 ) obtained from the study of Heller and Dobrzynski 
(ref. 8). The latter study presented far-field trailing edge noise spectra 
only as a function of frequency. However, the value of 6 “5 mm, presented 

in reference 8 to nondimensionalize the surface pressure spectra, can also 
be used to form a nondimens ional Strouhal number for the far-field radiation. 
Based on the spectrum peak occurring at f = 200 Hz when U = 50 m/sec (see 
fig. 14, ref. 8), the peak Strouhal number becomes St(6 ) = 0.02. This is 
approximately a factor of 4 to 5 smaller than the peak Strouhal value 
reported in reference 10. 

The Strouhal peak associated with the measurements of Heller and Dobrzynski 
(ref. 8) agrees with the average value determined from the present study. If 
St = 0.11, based on the average value of 6, then using 6 = 6/8, the Strouhal 

number based on displacment thickness is St(6*) “ 0.014. This is similar to 
the value of St (6 ) = 0.02 calculated for the study of Heller and Dobrzynski. 


Generalized Scaling Law 

Approach - The general conclusion from figures 29 to 31 is that for a 
given airfoil geometry, 1/3 octave band trailing edge noise spectra collapse 
onto approximately a single curve when normalized by the scaling law given in 
equation (66^ . It now remains to determine if the normalized spectra obtained 
from the different airfoil geometries can be represented by a single spectrum 
function, F 2 (St), and a constant, K^. For this purpose, the spectra in figures 
29 to 31 were plotted on the same graph as shown in figure 32. Only selected 
spectra are shown to avoid crowding the data. 

Development of Equation - The success of normalizing the spectra from 
different airfoil geometries demonstrates that there exists a universal 
trailing edge noise spectrum. To quantify the spectrum for noise predictions 
requires analytically describing the functional form for Fj^(St). This can be 
achieved by fitting an equation to the nondimens ional spectrum such as the solid 
solid line in figure 32. In this case, the analytical equation corresponds to 
the function presented by Fink in reference 6 and 7 for full-scale airframe 
trailing edge noise. The function is currently used in the NASA ANOPP Total 
Aircraft Analysis. The equation describing the relationship between 1/3 
octave band sound pressure levels and the overall sound pressure level is 
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(67) 


F,(f) = SPL,/3 -OASPL 



where is the frequency associated with the spectrum peak. Replacing 

frequencyoy Strouhal number gives the solid line in figure 32. In this case, 
^^MAX ~ selected based on the average peak Strouhal value determined 

from the different data sources described in the table in the previous sub- 
section. In addition, a value of St 


the data points in figure 32. 


MAX 


0.1 provided the best agreement with 


It now remains to determine the relationship between the spectrum shape 
and the absolute sound pressure level. The link between these parameters is 
controlled by the unknown constant, K-j^, in equation (66). Based on the decibel 
level corresponding to the value at the maximum of the solid line shown in 
figure 32, F^(St = 0.1) + = 84.0. But F^(St = 0.1) = -9.16 dB in equation 

(67). Thus, = 93.2 dB. 


The final equations for predicting both overall sound pressure levels and 
1/3 octave band sound pressure levels are then 


OASPL = 50 loQio + '0 'og,o ^ + 10 log,o D + K,, 

e 

K, = 93.2 , U in knots 


( 68 ) 


and 


SPU„, . 0»SPU * 10 l 09 „ |o6,3(^)" l^r* ° 


(69) 


^"max “ 

Effects of directivity have been included in equation (68) to generalize 
the results for application to all radiation angles. It should be emphasized 
that the constant, K^, was determined from various isolated airfoil data sets 
with the airfoils at a finite angle of attack (a % 6®). Thus, in the strictest 
sense, the above equations cannot be applied to the a = 0* test conditions. 
Fortunately, changes in angle of attack can be accounted for by introducing 
the correct value of 6, as shown by the good collapse of data in figure 21. 
Equation (68), therefore, is a general equation. 

Comparison With Existing Scaling Law - The above formulation can now be 
compared to the equations developed by Fink (refs. 6 and 7) to predict air- 
frame trailing edge noise. The functional dependence on M, 6, and St is 
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identical to the formulation given by Fink. However, the directivity function 
used by Fink for a wing with zero sweep-back angle was different from that 
given by equation (32) . Fink proposed that 


□ (^e* t) " ) 


(70) 


Here 6g is the angle between the flight path and the observer position measured 
from the approach direction. Recall that 9g in the present study is measured 
from the opposite direction so that cos^(6g/2) in the above equation becomes 
sin^(0g/2). Thus, equation (70) is similar to the numerator of equation (32) 
but lacks the Mach number dependent terms. This is because the directivity 
pattern given by equation (70) represents a low Mach number result which was 
based on the early theory of Ffowcs Williams and Hall (ref. 16). In contrast 
the present study provides a general expression given by equation (32) . This 
functional form included in the present scaling law, is considered to be the 
correct formulation for trailing edge noise. 

Another difference between the present scaling law and the result of 
Fink becomes evident when comparing the value of the constant Kjl* Equation (68) 
uses a value of ® 93.2 dB while Fink reported a value of = 101.5 dB. The 
latter value is based on measurements using microphones on posts 1.2 m above 
the ground. Accounting for ground reflections, which cause measured levels 
to be approximately 3 dB above free-field conditions at high frequencies, 
results in = 98.5. Hence, the value of associated with the scaling law 
in references 6 and 7 was high by approximately 5 dB when compared to the 
present study. This difference was evident earlier when Schlinker compared 
his trailing edge noise measurements with the prediction developed by Fink. 

The comparisons reported in reference 5 indicated approximately a 4 dB 
difference between measured and predicted values. 
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TRAILING EDGE NOISE PREDICTIONS 


Objective 


Two different trailing edge noise prediction methods were developed in 
the previous sections of this report consisting of a fundamental theoretical 
approach with no adjustable constants and a scaling law approach based on 
experimental data. Each approach modeled the noise as radiating from a 
stationary two-dimensional isolated airfoil segment. 

It now remains to assess the accuracy of the two prediction methods. 
Following this the stationary two-dimensional isolated airfoil formulation 
will be extended to the rotating coordinate system. Transformations de- 
scribed in the section titled Theoretical Formulation of the Trailing Edge 
Noise Mechanism will be used for this purpose. Sound pressure levels will 
be calculated and compared to specific full-scale helicopter test data. 


Assessment of Isolated Airfoil Noise Prediction Procedures 

Scaling Law Approach - It would appear that selecting values for 
StMAx based on experimental data removes any discrepancy between the 

scaling law prediction and experimental measurements. Supposedly this would 
provide an accurate scaling law method for predicting the trailing edge noise. 
The accuracy of the scaling law approach, however, is limited by how closely 
the initial data collapses onto a single curve. Although the solid line in 
figure 32 represents the average of the experimental data, the curve should 
be compared to the individual data sets used to generate this figure. 

Figure 29 compares the scaling law with the experimental results of the 

present study. Measurements at Strouhal numbers below unity show reasonable 

agreement with the solid line representing equation (70) . The Strouhal nvnnber 

associated with the spectrum peak is close to the value of St = 0.1 

MAX 

Data from the study of Schlinker (ref. 5) is compared to the scaling law 
in figure 30. Again the sound pressure levels are in good agreement with the 
measured values while the spectrum peak occurs near the value of Stj,j^ selected 
for the scaling law. Comparing the scaling law with the measurements of Brooks 
and Hodgson in figure 31 indicates a difference in Strouhal number associated 
with the spectrum peak. This difference was not apparent in figure 32 due to 
the heavy concentration of data points. 
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The detailed comparison available in figures 29 , 30 , and 31 emphasizes 
that the scaling law represents the average of the data points in figure 32 
and specific data sets can be different from the average line. Such differ- 
ences define the accuracy of the scaling law approach when compared to the 
isolated airfoil studies conducted by different Investigators. 

Theoretical Approach - Figure 33(a) compares the first principles trailing 
edge noise theory given by equations (21) and (22) with the measurements 
obtained from the present study. Note that the power spectral density form- 
ulation in equation 21 has been converted to one-third octave band form. Also 
the value of 6* needed in equation (22) was determined from figure 19 using ‘ 
the flat plate boundary layer approximation 6* * 6/8 * 6^/(1. 07) (8) . Recall 
that the accuracy of modeling the flow field as a flat plate boundary layer 
was discussed earlier in the section titled Airfoil Boundary Layer Character- 
istics. Additional discussions are given in Appendix D. 

A comparison of the predicted and measured spectra in figure 33(a) shows 
good agreement at high frequencies but poor agreement at low frequencies. 

A similar comparison with the narrowband spectra of Brooks and Hodgson, given 
in figure 34, also shows poor agreement. In the latter case the value of 6 
needed as an input to the present prediction was based on Table IV and the 
relationship 6 * = 6/8. 

Predictions in figure 33(a) and 34 are based on a flat plate surface 
pressure field model with the spectrum shape, S , given by equation (22). 

It was postulated that better agreement at low frequencies could be obtained 
by replacing the flat plate convecting surface pressure model with the pres- 
sure field existing at the trailing edge of an airfoil. Justification for 
this approach was based on the normalized surface pressure spectrum measure- 
ments of Yu (ref. 9) as well as Brooks and Hodgson (ref. 10). The experi- 
mental results obtained from reference 9 are plotted in figure 35(a) as a 
function of m for two freestream velocities. These measurements correspond 
to an airfoil at zero degrees angle of attack. Also shown in the same 
figure is the normalized surface pressure spectrum for a flat plate given by 
equation 22. 

Comparison of the airfoil and flat plate spectra in figure 35(a) indicates 
a definite difference in the spectrum shape and absolute level for these two 
geometries. Qualitatively, the fluctuating surface pressure spectrum near 
the airfoil trailing edge is seen to be about 7 dB higher than the zero pres- 
sure gradient flat plate result. This contrast is presented quantitatively 
in figure 35(b) where the abscissa represents the difference between the 
average of the airfoil data (see figure 35(a)) and the flat plate spectrum. 
Representing the difference function in figure 35(b) as DIFF (u) results in 
a modified expression for the surface pressure .spectrum given by 
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( 71 ) 


i ' (w.o) = (i -u f : 


Applying S' to the Isolated airfoil noise prediction provides the 
modified theoretical prediction shovm in figure 33^) for trailing edge noise 
measurements conducted during the present study. A comparison of the modified 
theory and experiment now shows better agreement at low frequencies. In 
particular, predicted noise levels are within 3 dB of measured levels over 
most of the spectrum range from M = 0.1 to 0.5. The need to include the air- 
foil surface pressure spectrum has, therefore, been demonstrated. 

A similar conclusion is obtained from figure 34 where the modified theory 
is compared to the measurements of Brooks and Hodgson. The improvement in 
the predicted sound pressure levels is evident when the circle and square symbols 
are compared at M = 0.2. 

A final demonstration of the ability to calculate trailing edge noise 
radiation using the modified first principles theory is given in figure 36. 

Here the present analysis is compared to the data reported by Schlinker 
(ref. 5). Again the comparison between theory and experiment shows good 
agreement in absolute sound pressure levels. 

Figure 34 also shows the trailing edge noise prediction calculated by 
Brooks and Hodgson for their oim data. A brief discussion of the difference 
between their prediction and the present calculation is warranted. The pre- 
diction in reference 10, which is based on the measured surface pressure 
spectra and characteristic length scales, does not require a value of 6 as 
an input. In contrast, the present approach uses an analytical model for the 
surface pressure characteristics and requires a knowledge of 6* which is 
calculated from flat plate boundary layer theory. As noted in Appendix D, 
the calculated 6* value differs from the measured value reported by Brooks 
and Hodgson by a factor of approximately 2.5. It is not clear which value 
should be used. However, Appendix D discusses the above described differ- 
ences in 6 and presents the reasons for using the value of 6* obtained 
from flat plate boundary layer theory. 

Modification of the first principles trailing edge noise theory focused 
on replacing the expression for Sqq on a flat plate with a measured airfoil 
surface pressure spectrum. A similar modification was considered for the 
spanwise length scale, £y, which is derived in equation (19) from the span- 
wise cross-spectrum Sqq (o),yo). Based on the comparisons of Sqq and S'qq the 
functional form of the parameter, iy, would not be expected to be the same 
for an airfoil surface and a flat plate. However, Brooks and Hodgson showed 
that. Indeed, the normalized airfoil and flat plate spectra are Identical. 
Ck)nsequently, if was not necessary to alter the expression for iy and the flat 
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plate boundary layer length scale was considered to be representative of 
the flow field modeled here. 


Comparison of Scaling Law and Theoretical Approach - It is worthwhile 
to compare the analytical and the scaling law approach for calculating trail- 
ing edge noise. For comparison purposes the theory will be simplified by re- 
moving the dependence on the ratio of acoustic wavelength to airfoil chord. 
Letting the ratio of chord-to-wave-length approach infinity, and letting 
= M in the expression for £ in equation (21) gives for the far-fleld 
narrow band sound pressure level 

Spp(x,o;)=^^ D(0e.<#*) Sqq(w.O) (72) 


where D is given by equation (32 ). Note that the modified surface 

pressure spectrum given by S’ in equation (71) has replaced the flat plate 
surface pressure spectrum terS^in equation (21). 


Putting this in third octave form using equation (24) gives approx- 
imately 

V II SB* 

SPL,„ • C 109,0 0(9...^)F(0)DIFF(B)J 

+ 126.8 0.1 < CO < 20 


where F(o)) represents the frequency dependence given by equation (22b). Here 
the value K follows strictly from the theoretical noise model proposed and 
the values of the parameter and obtained after assuming that surface 

pressures near the trailing edge can be modeled by equation (71) and (24). 

No adjustment of the value of K exists in equation (73). It should be noted 
that there is a factor of 8tt difference between equations (72) and (73). This 
change occurs because a factor of 2 is needed to account for both sides of 
the airfoil, a second factor of 271 accounts for changing from an w spectrum 
to a frequency spectrum, and an additional factor of 2 is needed to convert 
a two-sided spectrum ranging from o) = -<» to co = to a one-sided spectrum 
ranging from o) = 0 to co = -H». 


Before the above theoretical result can be compared to the scaling law 
equation it is necessary to convert the equation to a form similar to that 
given by equation (69) . This requires replacing the Mach number dependence by 
velocity and letting 6 =6/8. Equation (73) then becomes 

SPL,/3 = 50 log,Q ^lookf^ots) {r^) 10910(5) 

+ 10 IoQk) [f(0)DIFF (uJ)] + K (7A) 


O.K QJ < 20, K = 76.6 
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Now the constant, K, contains the conversion which changes Mach number 
dependence to a velocity dependence in knots. The constant also contains 
the factor of 8 which links 6* and 6. 


For comparison to the above expression, equations (68) and (69) are 
combined into one similar equation which is 


spl ,/5 = 50 log , 0 (|55^) + 10 log„ p + lo log 


10 


( 0 ) 


'0'°9lojo-613(s^j‘‘[(s5^)’^+0.5 


(75) 


+ K. 


where = 93.2. 

The form of equations (74) and (75) now permit direct comparison. The 
parametric dependence on U, 6, and D (6g,(j)) is similar. However, the peak 
value in addition to the frequency dependence appear to be different in each 
equation. These differences are best evaluated by isolating the frequency 
and amplitude dependent terms in a form similar to equation (66), thereby, 
permitting direct comparison of the normalized spectrum shape. The result 
is shown in figure 37 where the abscissa represents either: 

a) for the first principles theory, F(St)DIFF (St)+K 

b) for the scaling law, F 2 ^(St)+K^ 

Figure 37 shows that the spectrum peak amplitudes for the two prediction 
methods are approximately the same. However, the Strouhal numbers associated 
with each peak differ by a factor of two even though the spectrum shapes are 
similar (see shifted curve in figure 37). In the case of the scaling law the 
peak occurs at St = 0.1 while the theoretical result gives a value of approx- 
imately St = 0.2. This explains why the peak Strouhal number in the analytical 
prediction shown in figure 33(b) occurs near St = 0.2 while the scaling law 
(based on figure 32) peaks near St = 0.1. 

Siammary and Evaluation - The above discussions assessed the accuracy of 
the scaling law and theoretical approach for predicting trailing edge noise 
levels. Based on the comparisons between predictions and experiment both 
methods give approximately the same sound pressure level at the spectrum peak 
and predict approximately the same spectrum shape. However, the fundamental 
theory predicts a larger value for the Strouhal number associated with the 
spectrum peak. Thus, for the purpose of predicting trailing edge noise 
radiation from rotating blades the scaling law is considered presently to be 
the most reliable. This conclusion applies for a range of Mach numbers up 
to a maximum of M = 0.55. Furthermore, this result applies to Reynolds numbers 
representative of full-scale helicopter main rotor operating conditions. 
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For those cases in which the Mach number is above the highest value 
investigated here the fifth power Mach number scaling continues to be 
applicable. This follows from the fundamental theory which includes the 
compressibility effects occurring at high Mach numbers. Based on equation 
(72) SppaMU^. For a rotor operating in a constant atmospheric propagation 
speed given by c^, SppaM^. Thus, the Mach number dependence is the same 
at low speeds and high speeds. 

Trailing edge noise directivity dependence in the scaling law is also 
considered to remain the same at high rotor Mach numbers. This follows again 
from the theoretical directivity function which was incorporated directly into 
the scaling law. Similarly the spectrum shape is considered to remain unaltered 
based on the theoretical normalized spectrum. Recall that the analytical and 
empirical spectrum shapes given by F (St) 'DIFFCSt) and F^(St) have approximately 
the same shape although they peak at a different Strouhal number (see Fig. 37). 

There are two limitations to the applicability of the trailing edge noise 
prediction procedures developed here. The first requirement is that the air- 
foil section Mach number must be below the condition at which impulsive noise 
is generated due to locally supersonic flow. In this case the shock boundary 
layer interaction increases the trailing edge boundary layer thickness making 
the flat plate boundary layer calculation unreliable for determining the value 
of 6 needed as an input to the scaling law. It should be noted that based on 
the scaling law dependence trailing edge broadband noise would increase sig- 
nificantly under this operating condition. 

The second requirement applies to model noise tests and the location of 
the laminar to turbulent transition point on the airfoil. To simulate full 
scale trailing edge noise radiation it is required that turbulent boundary layer 
develop over a sufficient distance to ensure an equilibrium boundary layer at 
the airfoil trailing edge. Otherwise the surface pressure spectrum field con- 
vected over the trailing edge does not possess a universal spectrum shape and 
becomes a function of the turbulent flow distance on the airfoil surface. 

The above criterion is best satisfied by maintaining the same Reynolds number 
as that expected for a full scale geometry. The length scale in the Reynolds 
number calculation would be based on the distance from the transition point to 
the airfoil trailing edge. 

It should be noted that differences between the present theoretical 
approach and the experimental results should not preclude further development 
and application of the theory. The analytical model of surface pressure 
spectra convecting past the airfoil trailing edge is basically correct and 
can be shown to be similar to the formulations of other investigators. The 
source of the difference in spectrum peak location is, therefore, not the 
acoustic source model but rather the surface pressure spectra data needed as 
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an input to the prediction. The lack of agreement can be traced to 
insufficient details in the surface pressure field description. Specifically, 
the pressure spectrum, S’ , is presently based on measurements obtained for 
an airfoil at a=0^. Differences between the predicted and measured spectrum 
peak Strouhal number could possibly be removed by incorporating the influence 
of angle of attack on the surface pressure spectrum. Measurements of both 
Yu (ref. 9) as well as Brooks and Hodgson (ref. 10) have shown that increas- 
ing a increases the low frequency portion of the suction side surface pressure 
spectrum by approximately 3 dB. Conceptually this would increase the absolute 
sound pressure level of the low frequency side of the predicted spectrum in 
figure 37 sufficiently to shift the broad spectrum peak by the one octave in 
figure 37. This could result in total agreement of the first principles trail- 
ing edge noise theory with the experimental results. 

Accurate predictions of the trailing edge noise mechanism may, therefore, 
require using measured surface pressure data obtained near the airfoil trail- 
ing edge. Such detailed measurements were conducted by Brooks and Hodgson 
at low Mach numbers (M^O.2) to provide the necessary experimental input to 
their theoretical noise prediction. This approach should be extended to the 
higher Reynolds numbers associated with full-scale helicopter rotor operating 
conditions. If it can be demonstrated that the surface pressure data at higher 
Reynolds numbers continues to collapse onto a universal curve, then this curve 
could be used as a generalized input to the noise prediction. 


Rotating Blade Noise Predictions 

Approach - For the purpose of predicting trailing edge noise from rotating 
blades the scaling law was transformed to the rotating frame. The method 
employed assumed that the spectrum of a given blade segment is, at any partic- 
ular instant, given by equation (75). The local relative velocity of the rotor 
blade segment determined the boundary layer details and the resulting noise 
spectrum. The 1/3 octave band spectrum shape, given by equation (75), was 
converted to a narrowband power spectral density expression by accounting for 
the filter bandwidth. The spectrum was permitted to vary as the rotor moved 
about the azimuth. 

To obtain the final spectrum, the instantaneous spectrum was averaged 
around the azimuth, together with a weighting which accounted for retarded 
time effects. Details of the analysis were described earlier in the section 
titled Theoretical Formulation of the Trailing Edge Noise Mechanism. 

Input to Noise Prediction - Equation 75 requires for each spanwise segment 
of the rotor blade a knowledge of a) the local freestream velocity, b) the 
average of the suction side and pressure side boundary layer thickness, and 
c) the observer position relative to the rotor plane. With the exception of 
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the boundary layer thickness, 6, the above variables are easily determined. 

To avoid using complicated computational fluid dynamics methods for determin- 
ing 6 this parameter was estimated using the flat plate boundary layer cal- 
culations. This approach was justified earlier by a comparison of the flat 
plate calculation with the average of the pressure and suction side values 
of 6 for the rotor blade segment tested in the present study. The comparisons 
showed good agreement at small angles of attack (see fig. 19). 

Recall that the differences between flat plate boundary layer theory and 
experiment were less than 20 percent. Based on the trailing edge noise 
scaling law the estimated difference in noise radiation would be less than 
0.8 dB. This statement imples that helicopter rotor trailing edge noise 
predictions based on flat plate boundary layer characteristics are accurate 
within several decibels. Greater accuracy, and therefore, detailed boundary 
layer information, was not required since the objective of the study was to 
assess the relative importance of this noise mechanism compared to the other 
operative mechanisms on a full scale helicopter. 

Since the abrasion strip on full scale helicopter rotors is known to 
trip the surface boundary layer close to the blade leading edge the flat plate 
boundary layer calculation was initiated at the blade leading edge. This is 
expected to introduce only a small error in the calculation of 6. 


Test Case - The trailing edge noise prediction method was tested against 
data measured during a full-scale helicopter noise flyover. Details of the 
helicopter main rotor and tail rotor are tabulated below. 


Main Rotor Tail Rotor 

Number of Blades 
Rotor Radius 
Blade Chord 
Rotor Rotation Rate 


4 

6.7 m 
0.4 m 
307 rpm 


4 

1.2 m 
.17 m 
1689 rpm 


The main rotor blade geometry was identical to the rotor blade segment tested 
during the two-dimensional isolated airfoil investigation conducted in the 
present study. 

Data presented here correspond to a level flyover condition of 36 or 
72 m/sec and an altitude of 76 m. Acoustic spectra were processed when the 
helicopter was in approximately an overhead position so that 0^ = 90®. Data 
was obtained using microphones on posts 1.2 m above the ground. To account 
for the ground reflection effect 3 dB were added to the free field trailing 
edge noise prediction. Weather conditions during the helicopter flyover 
are tabulated below. 
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Wind speed - zero velocity (calm) 

Temperature - 77®F 

Relative Humidity - 90% 

The effect of atmospheric attenuation due to humidity was accounted for by 
applying a correction to the measured data. 

Figure 38 shows the narrowband (12.5 Hz bandwidth) total helicopter noise 
spectrum measured during the flyover. Corrections for ground reflections and 
atmospheric attenuation have not been applied here. The objective of pre- 
senting this figure is to demonstrate that the flyover signature is mainly 
broadband noise above 1 KHz. Apparent discrete tone peaks in the spectrum 
are due only to the short averaging time (1/2 sec) used during the spectrum 
analysis. Although a longer averaging time would have provided a smoother 
spectrum, spatial directivity information would then be averaged. For example, 
the 1/2 second averaging time employed in the above described measurements 
resulted in averaging angular information within - 11® about the 90® overhead 
position during the 72 m/sec flyover. Fortunately, directivity information 
varied slowly at this angle as shown by the three different spectra super- 
imposed on figure 38. Each spectrum represents the output from a 1/2 sec 
spectrum analysis obtained during the flyover. One spectrum corresponds to 
the 90® overhead position while the other two spectra represent the preceed- 
ing and following spectra. 

Figure 39(a) shows the same flyover condition measured with a narrowband 
spectrum in figure 38 converted to 1/3 octave band analyzer (circle symbol). 
The square symbol in figure 39(a) represents the same data with the atmospheric 
attenuation correction added to the measurement. Additional corrections to 
account for the Doppler shift and retarded source position were not necessary 
in the overhead position. 

Also shown in figure 39(a) is the scaling law prediction for the helicop- 
ter main rotor trailing edge noise. Comparison of the total helicopter noise 
spectrum and the trailing edge noise spectrum demonstrates that at high fre- 
quencies trailing edge noise can be a significant noise mechanism. In fact, 
at some frequencies predicted and measured noise levels are almost identical. 
Below 2 KHz other broadband noise mechanisms begin to dominate over the trail- 
ing edge noise as shown by the divergence of the predicted and measured curves. 

Figure 39(b) shows a second measured spectrum at a different flyover 
speed. Again the comparison between predicted 1/3 octave band levels and 
measured spectra is good. 

Also presented in figure 39 is the first principles prediction for the 
main rotor tailing edge noise. The spectrum peak amplitude is close to the 
measured levels although the frequency associated with the peak occurs at a 
value about the dominant portion of the measured broadband noise spectrum. 
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This latter feature was expected from the comparison of the scaling law 
and the first principles trailing edge noise predictions in figure 37. It 
is felt, however, that further development of the airfoil surface pressure 
spectrum model will provide better agreement with measured data. 

A comparison of the predicted trailing edge noise for the helicopter 
main rotor and the measured total helicopter noise spectrum has been pre- 
sented in figure 39. Although noise from the tail rotor was also present 
this contribution was expected to be weak. This is due to the smaller bound- 
ary layer thickness and rotor span values which are linear inputs to the 
analysis according to equation (75). Also, the smaller boundary layer 
thickness on the tail rotor causes the spectrum peak to shift to a higher 
Strouhal number. 

Summary and Evaluation - The importance of the helicopter trailing edge 
noise mechanism has been verified by direct calculation of the associated 
noise levels during an aircraft flyover. Predicted noise levels close to the 
total helicopter noise spectrum demonstrate that trailing edge noise must be 
accounted for in future noise prediction procedures. 

Since trailing edge noise is controlled by the inherent characteristics 
of the rotor blade boundary layer, it is difficult to eliminate this noise 
mechanism. Delaying the boundary layer transition to reduce 6 in equation 
(75) is a possibility but the benefits are marginal due to the linear depen- 
dency on boundary layer thickness. Efforts to shift the blade surface pres- 
sure spectrum to a higher frequency to increase atmospheric attenuation re- 
presents a possible but difficult task of controlling details of the turbulent 
boundary layer pressure spectrum. 

The reduction of helicopter rotor trailing edge noise is, therefore, 
limited to changes in the aerodynamic operating conditions. As suggested 
by the prediction methods presented here, a. reduction of tip speed would 
provide the greatest improvement. Reduced speeds would, presumably, re- 
quire increasing the number of rotor blades to maintain a constant rotor 
thrust. The resulting increase in trailing edge noise would, however, be 
less than the decrease achieved by changing the tip speed. 

With the exception of the above described noise reductions, trailing 
edge noise represents the limiting acoustic source mechanism for the total 
helicopter noise spectrum. This conclusion is similar to the conclusions 
obtained in previous studies of the noise generated by fixed wing aircraft. 
This noise, referred to as "airframe noise", has set the noise limit for 
such aircraft. Based on the results of the present study, a similar noise 
limit exists for helicopter rotor broadband noise. 
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Reductions in high frequency helicopter broadband noise have also been 
achieved by changing the rotor blade tip shape. The source mechanism, in 
this case, is not the attached boundary layer trailing edge noise investigated 
in the present study. Instead, the noise is generated by convection of the 
three dimensional tip vortex flow field over the blade trailing edge. George 
(ref. 50) provided an approximate model of the noise generated by this 
separated and highly turbulent flow field as it convects over an airfoil trail- 
ing edge. His predictions indicated that blade tip vortex noise can be 
comparable to trailing edge noise. This result explains the broadband 
noise reduction observed by some investigators after changing the airfoil 
tip shape. 


\ 
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CONCLUSIONS AND RECOMMENDATIONS 


A. Full-'Scale Rotor Noise 

A.l. Extension of a validated isolated airfoil trailing edge noise scaling 
law to the rotating blade case demonstrated that trailing edge noise from a 
full-scale helicopter contributes signif icantly to the total broadband noise 
spectrum at high frequencies. This noise mechanism is expected to control 
the minimum rotor noise level. 

A. 2. The helicopter rotor trailing edge noise scaling law for subsonic tip 
speeds was validated using two-dimensional isolated airfoil results at full- 
scale Reynolds numbers. Knowledge of the helicopter altitude, speed, angular 
position, rotor tip speed, blade number, rotor chord and span, and rotor 
boundary layer thickness are the only parameters required in the analysis. 

Based on two-dimensional isolated airfoil results, helicopter rotor boundary 
layer thickness can be estimated using a flat plate turbulent boundary 
layer calculation. This simplifies the flow field details needed as an 
input to the prediction procedure. 

A. 3. A first principle theory for trailing edge noise predicts the correct 
absolute sound pressure level but does not predict the Strouhal value identi- 
fying the spectrum peak. The basic noise model of a frozen pressure pattern 
convecting past a trailing edge, however, agrees with the acoustic source 
models developed by other investigators. Differences between theory and 
experiment are believed to be due to the airfoil fluctuating surface pressure 
field model used as an input to the analysis. Accurate predictions require 
using measured surface pressure data obtained near the airfoil trailing edge. 

A. 4 Future efforts should be directed toward detailed airfoil surface pressure 
measurements at Reynolds number representative of full scale helicopter opera- 
ting conditions. If it can be demonstrated that the surface pressure spectra 
data at high Reynolds numbers collapse onto a single nondimens ional curve, 
then this curve can be used as an input to the noise prediction. 


B. Trailing Edge Noise From a Local Blade Segment 

B.l. Based on two-dimensional airfoil results, trailing edge noise radiation 
patterns from a local blade segment are well predicted by the theory developed 
in the present study. Directivity patterns are independent of angle of attack 
and acoustic source Strouhal number. 


64 



B.2. Noise dependence on local section Mach number varies as M^ and is 
independent of angle of attack and acoustic source Strouhal number. 

B.3. Local rotor blade angle of attack influences only the low Strouhal 
number portion of the spectrum. Increasing angle of attack increases sound 
pressure level by only a few decibels near the spectrum peak. The observed 
trend is explained by the increase in boundary layer thickness on the suction 
side of the airfoil. 

B.4. Trailing edge noise 1/3 octave band spectra from a local blade segment 
can be approximated by a universal spectrum shape. A scaling law capable of 
predicting the absolute spectrum levels depends only on Mach number, boundary 
layer thickness, and observer location. 

B.5. Boundary layer thickness input data to the scaling law can be estimated 
using flat plate turbulent boundary layer calculations. Such calculations 
approximate the average value of the pressure side and suction side boundary 
layers. The average value controls the peak Strouhal number of the universal 
spectrum curve. 
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APPENDIX A 


Directional Microphone System 


Directional Sensitivity and Spatial Resolution - Although the directional 
microphone system is most sensitive to sources located on its aiming axis, it 
does sense sources at off axis directions defined by the angle 4)^ in figure 12. 
The response of the system to such off axis sources is controlled by diffrac- 
tion at the circular aperture of the reflector. For the present geometry 
the diffraction phenomenon is governed by the nondimensional parameter 
n = (irDf sin(()o) /C q where D is the reflector aperture, f is the acoustic fre- 
quency, and Cq is the sound speed. Figure 40 shows the measured normalized 
response from the focal point microphone plotted as a function of n* The 
response function, H, decreases as 4)^ increases demonstrating the decrease 
in sensitivity to off-axis angles. This response, often referred to as the 
spatial discrimination capability, has been shown to agree with the Fraunhofer 
diffraction pattern for a circular aperture system (see Fig. 40 ). Such 
agreement permitted using an analytical expression for the response function 
in the data reduction. 

Figure 40 demonstrates that diffraction causes the system to have a finite 
spatial resolution instead of the ideal delta function resolution. Thus, 
while the system is most sensitive to sources located on its aiming axis, it 
does sense sources at off axis directions. The "sharpness” of the resolution, 
however, increases with the acoustic frequency, f, which appears in the numer- 
ator of the parameter, n- Thus, for a fixed off-axis source position defined 
by 4 ^q, the response H decreases as f increases. 

One parameter which quantifies the directional microphone frequency 
dependent spatial discrimination is the diffraction pattern half-width, AW. 

The half-width, defined as the spatial distance between the negative 3 dB 
points in the diffraction pattern of figure 40 corresponds to the points 
r) = - 1.66. A 3 dB decrease in sensitivity to off-axis sources, therefore, 
occurs for displacements given by x = 1.66 c^R/fD where x corresponds to the 
off-axis distance shown in figure 12. Here sin x/R has been employed in 

the expression for n. The resolving half-width is then AW = 2x = 3.32 c^R/fD. 
This result is shown in figure 41 in addition to the experimentally determined 
value for a source situated at a distance of R = 2.81 m from the reflector 
The good agreement between theory and experiment confirms the accuracy of the 
Fraunhofer diffraction pattern description for the reflector spatial discrimina- 
tion characteristics. 
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Gain - The directional microphone system also has a significant gain. A 
gain exists because the large collecting area of the reflector focuses the 
sound onto the small focal point microphone. In an ideal system without dif- 
fraction the gain, G, depends on the solid angle subtended by the reflector 
and the focal point microphone so that G = (D/d)^ where D = 1.067 m and 
d = 0.635 cm. However, the frequency dependent diffraction characteristics 
smear out the focal point image so that only part of the acoustic energy 
incident on the reflector is concentrated on the focal point microphone. 

The directional microphone gain in, this case, can be determined 
experimentally (ref. 5). The sound pressure level generated by a point 
source of sound is measured with the reflector and with an omnidirectional 
microphone. Source- to-ref lector and source-to-omnidirectional microphone 
distances are identical during the measurements. The ratio of these two 
measurements represents the gain, G, 


G (f,R) = 


(f.R) 


(A.l) 


2 2 
where defines the focal point microphone sound pressure level and p^ 

represents the omnidirectional microphone sound pressure level. The gain 

measurement, which was repeated in the present study as a check on the 

directional microphone system, is shown in figure 41 for two source distances 

corresponding to R = 2.81 m and 2.07 m. The large gain was beneficial at 

high frequencies because turbulent boundary layer noise decreases rapidly 

with increasing frequency. 


The directional microphone gain must be determined before absolute 
source levels can be calculated. At low frequencies the gain function can 
be calculated explicitly (ref. 5) based on the Fraunhofer diffraction theory 
giving the result 


G(f,R) = ^ 


jTOi 

2RCq 


2R-Rc1 2 2 


(A. 2) 


which is plotted in figure 41. 
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Differences between the measured and theoretical gains at high frequencies 
in figure A1 are attributed to spherical aberrations at the focal point 
microphone. These characteristics are not included in the analytical calcula- 
tions of the gain. Thus, the experimentally determined gain was used at high 
frequencies for all data reduction. 
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APPENDIX B 


Forward Flight Effects 


Shear Layer Refraction - As shown in figure 13, refraction changes both 
the acoustic ray path as well as the path length of the sound waves arriving 
at a microphone station outside the flow. In addition, transmission of sound 
across the open jet shear layer changes Che divergence rate of the acoustic 
ray tubes. These latter changes in path length and divergence rate require 
correcting the measured sound pressure levels. 

The detailed analytical study and experimental verification by Schlinker 
and Amiet (ref. 47) formed the basis of the corrections applied to the present 
data. The corrections were used to convert measurements at the observer 
station (0), in figure 13, to the constant sideline station (D) on the non- 
refracted path. The corrected sound pressure level represents an equivalent 
measurement in the absence of the shear layer. In other words, the shear 
layer is considered to be located at infinity with uniform flow existing 
between the source and the far field microphone at station D. As concluded 
in the study of Schlinker and Amiet, the finite thickness shear layer can be 
replaced by an idealized vortex sheet. Specific details of the shear layer 
mean velocity profile are, therefore, not needed. 

An example of the magnitude of the angle and amplitude changes is given 
in figure 42 for an open jet Mach number of M = 0.5, corresponding to one of 
the operating conditions tested in the present study. The distance to the 
sideline microphone measurement station is given by h/y^^ = 0.14 where h 
represents the source-to-shear layer distance and is the sideline distance 
in figure 13. As shown by figure 42 , sound initially propagating at 
6^ = 80° inside the open jet is refracted to 0^. = 113° outside the airstream. 
Also, the sound pressure level measured on the resulting refracted path at 
station 0 in figure 13 should be increased by 4 dB according to figure 42. 
Note that the amplitude correction is always added to the measured data. 

Apparent Source Location - The above described refraction angle changes 
control the propagation of a single ray through the shear layer and arriving 
at the far field microphone station (0) in figure 13. In contrast, the 
directional microphone system senses refracted acoustic rays propagating 
within the angle, d0^, as shown in figure 14(a). Transmitted rays, A and B, 
each of which are refracted differently, define the angle, d0^, in the hori- 
zontal plane. Since the directional microphone system senses only the 
refracted rays outside the flow, it must be focused on the **apparent” source 
location defined by rays A and B. The intersection of these rays inside the 
flow controls this apparent source location in figure 14(a). 
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Position was determined analytically using the refraction angle 
correction controlling rays A and B. (See Appendix C for a detailed evalua- 
tion.) For each Mach number, M, and directivity angle, 0^, a unique 
apparent source position was calculated. As shown by the analytical deriva- 
tion, the apparent source location is independent of the distance to the re- 
flector surface. 

Once the apparent source location was calculated, the sideline distance, 
y^, for the directional microphone traverse could be selected. This distance 
was controlled by the apparent source-to-ref lector distance, R, in figure 12. 
Since detailed calibrations of the directional microphone spatial resolution 
and gain were only available for two specific values of R, the value of 
was selected to provide a distance of R - 2.07 m or 2.81 m between the 
reflector and the apparent source position. Recall that the value of R con- 
trolled the focal point microphone location r^, in figure 12. 

As shown by figure 14(a), the directional microphone senses sound within 
a solid angle. In the horizontal plane, rays A and B (fig. 14(a)) describe 
the limits of this solid angle in addition to determining the apparent source 
position. In the vertical plane (fig. 14(b)) rays C and D define similar 
limits of the solid angle. According to the arguments previously presented, 
the directional microphone system must be focused on the apparent source 
position defined by the intersection of rays C and D inside the flow. These 
rays experience identical refraction angle changes in the vertical plane 
in contrast to the horizontal plane where rays A and B are refracted differ- 
ently as they propagate through the shear layer. Thus, tracing rays C and D 
back into the flowfield resulted in a different apparent source position in 
the vertical plane when compared to the horizontal plane. Details of the 
analysis used to predict the apparent source position, Sy (see fig. 14(b)) 
in the vertical plane are included in Appendix C. 

Both S„ and S,, are situated on the reflector centerline but at different 
n V 

distances from the open jet shear layer. Figure 14(a) shows that is 
situated at a distance, a, from the shear layer while figure 14(b) indicates 
that Sy is located at a distance of b^* The significance of these different 
apparent source locations is best demonstrated by a numerical example. For 
M 0.5 and 0^ = 50®, the direct distance between and is 16 cm. 

The ability to sense the acoustic radiation from these different source 
regions depends on the directional microphone depth of field. This parameter, 
shown in figure 43 represents the frequency dependent decrease in the 
measured directional microphone output as a source is moved closer or further 
from the reflector. In each case the focal point microphone distance, r^, is 
fixed by the distance R shown in figure 12. Consequently, source position 
variations about R, represented by R ! A, result in a decrease in the measured 
response. 
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The different source locations given by and S_, represent the extremes 

V 

of the apparent source position in the open jet. Ideally, the directional 
microphone should have a large depth of field to sense the distributed 
apparent source region. Oddly, this contradicts the usual requirement for 
sharp spatial resolution demanded of directional microphone systems. 

The final criteria selected for focusing the directional microphone sys- 
tem on the apparent trailing edge noise source region was based on the distance 
to the midpoint between and S^. As shown in figure 44, the sideline dis- 
tance, yj^, was selected to provide *= R* Here R = 2.81 m or 2.07 m 

corresponding to the two distances for which detailed depth of field calibra- 
tions were available. Most important of all, with the system focused on 
the decrease in the system output due to the source at - AR/2, 

where AR represents the distance between and Sy, could not exceed 0.5 dB. 
This tight tolerance assured that the apparent source region was definitely 
within the reflector depth of field. All data reported in the present study 
satisfied this condition. 

It should be noted that as the measurement angle, 6^, approached the open 
jet axis, the distance between Sjj and Sy increased to the point where the 
depth of field was insufficient for sensing the entire source region at high 
Mach numbers. This situation was considered unacceptable and measurements at 
such conditions were not conducted at high Mach numbers for angles close to the 
open jet axis. Other investigators should be aware of this limitation when 
using a directional microphone to investigate acoustic source distributions 
at high Mach numbers in an open jet test section. 

Scattering Effect on Directional Microphone Measurements - The directional 
microphone spatial resolution and gain characteristics shown in figures 40 and 
41 apply only in the absence of flow between the acoustic source and the direc- 
tional microphone system. In reality, the directional microphone, situated 
outside the flow, collects the sound transmitted through the turbulent shear 
layer. Turbulent eddies in the shear layer can spatially scatter the sound 
sufficiently to reduce the measured acoustic source strength and broaden the 
source distribution pattern. For example, if the directional microphone 
were scanned past an ideal point source of sound inside the open jet the 
"apparent** gain of the system would be reduced while the diffraction pattern 
would be broadened. These changes must be accounted for since the data 
reduction procedure includes the directional microphone gain and spatial 
resolution characteristics. A summary of the effects of scattering on the 
measured sound pressure levels is provided here. 

To help understand the scattering effect on the measured directional 
microphone sound pressure levels, the following tentative physical explana- 
tion is presented. The argument is based on a point source of sound situated 
in the open jet potential core. Initially, a wavefront originating from the 
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acoustic source has a smooth spherical shape as shown in figure 45. While 
passing through the shear layer local scattering by the turbulence structure 
reorients small segments of the acoustic wavefront* The angle between the 
nomal to each wavefront segment and the reflector axis is no longer zero. 

Thus, the wavefront segment arriving at the reflector appears to originate 
from an off-axis source* Since the response to off-axis sources is diminished, 
the directional microphone output is decreased* In this manner, the effective 
gain of the system is reduced* At the same time, the diffraction pattern 
for the point source is broadened giving the appearance of a distributed 
acoustic source region instead of a point source inside the open jet* This • 
is similar to the appearance of a light observed through frosted glass. 

The size of the light appears to be much larger* 

The effects can be demonstrated quantitatively using the source distri- 
bution measurements shown in figure 46 * The experiment was conducted during 
the present study with the aid of a point source of sound located inside the 
open jet test section* In each case the directional microphone was scanned 
past the source providing 1/3 octave band source distribution curves like 
those shown for 10 kHz and 40 kHz* The abscissa in each figure represents 
the ratio of the directional microphone sound pressure level to the omni- 
directional sound pressure level. The latter microphone was placed at the 
directional microphone station corresponding to the peak in the measured 
acoustic source distribution* Apparent source- to-ref lector and apparent 
souce-to-omnidirectional microphone distances were identical. The ratio 
of the sound pressure level measurements at the source distribution peak 
represented the directional microphone gain defined earlier. 

At M = 0, scanning the directional microphone past the acoustic source 
provided the expected point source diffraction pattern and system gain* This 
is verified by comparing the measured half-width, AW in figure 46, with the 
theoretical value in figure 41 at the corresponding frequency* Also, a 
comparison of the gain represented by the peak values in figure 46 shows 
agreement with the previously measured values described in figure 41* 

As the open jet Mach number increased the peak amplitudes in figure 46 
decreased while the source distribution became broader. These changes were 
more pronounced at the higher frequency of 40 kHz since the acoustic wave- 
length was then comparable to the turbulent eddy size. Turbulence scatter- 
ing of sound has been shown to be stronger under this condition (ref. 47). 

The above described changes, represented by the parameters AMP and AW, 
were monitored over the range of 1/3 octave band center frequencies from 
1*25 kHz to 50 kHz. Measurements were conducted using the point source of 
sound at each Mach number and directivity angle investigated in the airfoil 
trailing edge noise study. As an example, the final effective gain, Gg, 
applied to the reflector measurements obtained at M = 0*5 and 0^ = 90° is 
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given by figure 47. Here the value of at f = 10 kHz, M = 0.5, and 
R = 2.07 m, is obtained from the corresponding peak value in figure 46. 

Also shown is the diffraction pattern half-width function, AW, which increases 
as frequency increases. These measured values obtained in the presence of 
flow, can be compared with the M = 0 measurements in figure 41. The difference 
in the directional microphone response demonstrates why scattering effects 
must be accounted for. Knowing the effective gain, G^, trailing edge noise 
sound pressure levels can be calculated from the directional microphone 
measurements. 
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APPENDIX C 


Apparent Source Position Beneath Shear Layer 

Because of the angle change of a sound ray as it passes through the 
shear layer, the apparent position of a sound source, from an acoustical 
standpoint, is not the same as its true source position. The apparent source 
position can be readily calculated. 

Consider rays in the plane defined by the vector normal to the two-dimen- 
sional shear layer and the flow vector and define this to be the xy plane (see 
fig, 48(a)). The actual source position is S and the apparent source position 
is The following geometrical relationships follow immediately: 



s = h/sin 6q 

(C.l) 


di = sd^c/sin ©c 

(C.2) 

Thus, 

d£ H hd 6^ /sin^ 6^ 

(C.3) 

Likewise, in terms of 

6^. rather than 0 
t c 



di = od0|/sin^ 9^ 

(C.4) 

Combining equations (C.3) and (C.4) 



sin^e, d6. 
a = h — -7^ 
sin^ 6^ dfl| 

(C.5) 

Equation (1) of Amiet 

(ref. 34) is 



tan dfs = C,/(/9^cos 6^ + M) 

(C.6) 

Taking the derivative 

gives 



. sinSt 

de, ■ ^3 ® 

(C.7) 
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which combined with equation (C.5) gives finally 


Q 



(C.8) 


For 6^ = 90®, the actual and apparent source positions will coincide, but at 
other angles they will differ. 

An apparent source position can also be calculated based on the ray paths 
in other planes. In particular, consider a plane perpendicular to the above 
plane and containing a ray in the xy plane after refraction by the shear layer. 
Thus, the plane makes an angle 0^ with the shear layer. This plane is shown 
in figure 48(b) which is basically the same as figure 11 of reference 53. The 
actual source position is not shown, but is again taken to be a distance h 
inside the shear layer. 

If dzj^ is the ray tube width just at the shear layer and dz 2 is the width 
at a distance of y^^ from the source, then equation (A. 19) of reference 53 gives 

^ = I + (y - l) 8. (C.9) 


From the geometry of figures 14(b) and 48(b) 

t| t| + (y, - h) CSC 
d2| " dZg 

Substitution for dz 2 /dZj^ from equation (C.9) gives 


(C.IO) 


t, = h/£, 


(C.ll) 


However, t^ is not the distance measured normal to the shear layer (fig. (14b)) 
since the plane considered makes an angle 6^ with the shear layer. The dis- 
tance normal to the shear layer is 


b| = t, sin 9^ - h 


sin0| 


(C.12) 


Again, for - 90®, = h, but in general bj^ h. Note that for 0^. ^ 90°, 
from equations (C.8) and (C.12), a b^ so that the apparent source position 
depends on which plane is being considered. 
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APPENDIX D 


Boundary Layer Thickness 

Boundary Layer Thickness Calculated for the Study of Schlinker (ref, 5) - 
For the NACA 0012 symmetric airfoil (22,9 cm chord) measurements reported by 
Schlinker, the boundary layer thickness was estimated using the flat plate 
boundary layer theory. It should be noted that a boundary layer trip was 
needed only on the pressure surface of the airfoil to eliminate discrete vor- 
tex shedding tones when the airfoil was at angle of attack. The trip was 
located at the 30 percent chord station. In contrast, natural transition was 
permitted to occur on the suction side of the airfoil. Based on the experi- 
mental data transition on the suction surface occurs at x/c = 0.05 when the 
airfoil operates at a = 6° angle of attack. A knowledge of the transition 
location on each surface was then used in the flat plate prediction for 6. 

The cumulative length over which the boundary layer developed corresponded to 
the distance from the individual transition points to the airfoil trailing 
edge. Contributions from the upstream laminar boundary layer were assumed 
negligible. The final value of 6 listed in Table III represents the average 
of the pressure and suction side boundary layers estimated for the study in 
reference 5. The credibility of such average value predictions was demon- 
strated by the close agreement between flat plate predictions and measurements 
obtained in the present study for an airfoil at small angles of attack (fig. 
19). Included in Table III are the airfoil geometric parameters. 

A check on the ability to predict 6 using the flat plate boundary layer 
theory is obtained by comparison with the NACA 0012 measurements reported by 
von Doenhoff (ref. 51). Logarithmic plots of the boundary layer velocity pro- 
files were obtained for a = 0° with natural transition occurring in the 
boundary layer. Considering the Re = 2.67 x 10^ case shown in figure 19 of 
reference 50, skin friction measurements indicated that transition occurred 
at approximately 50 percent chord when a = 0‘^. Calculating the equivalent 
flat plate turbulent boundary layer which starts developing at this station 
gives 6/c = 0.022 at the airfoil trailing edge. In comparison, the measure- 
ments of von Doenhoff indicate that the velocity profile at 97 percent chord 
asymptotically approaches unity near y/c = 0.02. Thus, the calculated flat 
plate result is in good agreement with the experimental result. In addition, 
reference 50 shows that at the 52 percent chord station 6/c = 0,0025 verifying 
that the laminar boundary layer upstream of the transition point represented 
only 10 percent of the final trailing edge boundary layer thickness. 

Boundary Layer Thickness Calculated for the Study of Brooks and Hodgson 
(ref. 10) - The trailing edge turbulent boundary layer thickness for the 
NACA 0012 airfoil (60.96 cm) chord was calculated using the flat plate boundary 
layer theory. The approach was similar to that used in calculating 6/c for 
the study of Schlinker (ref, 5) . Since a boundary layer trip was applied at 



the 15 percent chord station in the study by Brooks and Hodgson, the calculation 
was initiated at this station. 

Table IV gives the calculated value of 6/c in addition to the parameters 
describing the airfoil. The value of 6/c = 0.0166 represents the average of 
the pressure side and suction side boundary layer thickness values. 

It is worthwhile to compare the measured values of displacement thickness, 
6*/c, reported by Brooks and Hodgson with the value of 6/c calculated here. 

This comparison will be needed in a future discussion of the Strouhal value 
associated with the spectrum peak reported in reference 10. The present discus- 
sion will apply to the 6*/c data reported for the a = 0® angle-of-attack 
geometry. Data for the a = 5® angle of attack was not provided by Brooks and 
Hodgson. 

For the range of Mach numbers investigated at a = 0°, 6 /c was approxi- 
mately constant and equal to 0.0065 in reference 10. Then 6/6* 2! 2 based on 
the value of 6/c calculated above. This ratio is different from the typical 
value of 6/6 =8 measured in a two-dimensional flat plate boundary layer 

theory. One possible explanation could be the basic difference between the 
flat plate geometry used to calculate 6/c and the airfoil geometry for which 
6*/c was measured. This, however, would imply that the flat plate boundary 
layer calculation for 6/c is in error by a factor of approximately four. Yet, 
measurements obtained by von Doenhoff (ref. 52) showed good agreement with the 
flat plate calculation for 6/c at cx = 0° as noted in the previous subsection. 
Thus, it is possible that the value of 6 /c is high in the study reported by 
Brooks and Hodgson. 

A quantitative justification for the above statement is available from the 
recent study of D’Ambra and Damongeot (ref. 53) using a NACA 0012 airfoil. 

Their investigation of the airfoil fluctuating surface pressures and broadband 
noise presented frequency information in terms of a Strouhal number based on 
6*/c. Their displacement thickness was determined through sophisticated calcu- 
lations based on the potential flow field method of Garabedian and Korn com- 
bined with a viscous boundary layer calculation. The following simple relation- 
ship was obtained to evaluate the effects of the various parameters: 

Iog,o8* = -3.775+ \.\^ +0.12 M +0.8 log,QC + 0.775 a 

Here, x/c represents the distance from the airfoil leading edge to the station 
at which 6* is to be calculated. The airfoil chord length is specified in 
meters while a is given in degrees. 

The above equation is specific to the NACA 0012 airfoil and includes the 
location of the transition point and several other parameters. The equation 
represents the best curve fit derived from the theoretical computations and is 
applicable for the compressible flow range of 0.3 ^ M ^ 0.84 and -4® < a < 6®. 
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Extending equation 62 to the incompressible flow test condition given by M = 
0,2 in the study by Brooks and Hodgson gives 6*/c = 0.00247 for a = 0° . In 
comparison the value reported in reference 10 is 6*/c = 0.0065. Differences 
between these values suggest that the displacement thickness in reference 10 
may be large by a factor of approximately 2.5. 
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Table I Mach Number vs Propagation Angle, and Airfoil Angle of Attack 

GEOMETRIC 
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Table II Angle and Amplitude Changes Associated with Refraction and 
Retarded Source Position Corrections 


PROPAGATION ANGLE, e^, deg 



KEY TO PARAMETERS 
NUMBERS REFER 
TO LISTING IN 
EACH BOX. 


1) degrees 

2) degrees 

3) REFRACTION AMPLITUDE CORRECTION, dB 

4) AMPLITUDE CORRECTION FOR CONVERSION 
TO RETARDED SOURCE POSITION 
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Table III Average Boundary Layer Thickness Calculated for the Data of Schlinker 


VELOCITY 

M 

die 

m/sec 

30.5 

0.09 

0.0234 

45.7 

0.13 

0.0218 

61.0 

0.18 

0.0200 

76.2 

0.22 

0.0197 

91.5 

0.26 

0.0189 

106.7 

0.31 

0.0184 


AIRFOIL: NACA 0012 
CHORD: c = 22.9 cm 
SPAN: s = 53.3 cm 

SOURCE TO MICROPHONE DISTANCE: rg = 2.81 


Table IV Average Boundary Layer Thickness Calculated for the Data of 
Brooks and Hodgson 


VELOCITY 

m/sec 

M 

6/c 

38.6 

0.1 

0.0187 

69.5 

0.2 

0.0166 


AIRFOIL: NACA 0012 
CHORD: c = 60.96 cm 
SPAN: s = 46 cm 

SOURCE TO MICROPHONE DISTANCE: = 1.22 m 
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Figure 1 Simulation of Rotating Blade Trailing Edge Noise Radiation by a 
Two-Dimensional Airfoil Section 


87 










DIRECTIONAL 

microphone 

REFLECTOR 



91 


Figure 5 Open Jet Test Section Viewed from Far-Field Microphone Location 
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Figure 6 Section of Heiicopter Main Rotor Blade used for Isolated Airfoil Study 
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Figure 7 Rotor Blade Cross-Section 
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Figure 8 Location of Minimum Static Pressure Point on Airfoii Suction Side 




Figure 9 Facility Background Noise Sources 
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LENS EQUATION: 



Figure 12 Directional Microphone Operating Principles 
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Figure 13 Geometry for Transmission of Sound Through a Shear Layer 
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Figure 14 Apparent Source Position Effects on Directionai Microphone Alignment 
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Figure 16 Normalized Mean Velocity and Turbulence Intensity Profiles for 
Rotor Blade Section, a = *0.4° 
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Figure 16 — Continued 
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Figure 16 — Concluded 
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Figure 17 Normalized Boundary Layer Interface Distance Versus 

Mach Number and Angle of Attack for Rotor Blade Section 
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Figure 18 Normalized Mean Velocity and Turbulence Intensity Profiles for 
Rotor Blade Section, 6=7.B^ 
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Figure 18 — Continued 
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Figure 18 — Concluded 
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AVERAGE OF dJc OR 6/C 



Figure 19 Average of Normalized Boundary Layer Interface Distances on 
Pressure and Suction Side of Rotor Blade Section 
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Figure 20 Measured Axial Acoustic Source Distributions for Rotor Blade Section 
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Figure 21 Effect of Angle of Attack on 1/3 Octave Band Trailing Edge Noise 
Spectra for Rotor Blade Section at Various Mach Numbers, 6^-7 
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Figure 21 — Concluded 
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Figure 22 Comparison of Measured and Theoreticai Source Directivity as a 
Function of Mach Number for Rotor Biade Section, ct -*0.4^ 
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Figure 23 Comparison of Measured and Theoreticai Source Directivity for 
Rotor Blade Section, M =0.3, a = 7.6® 
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Figure 24 Variation of 1/3 Octave Band Traiiing Edge Noise Spectra with Mach 
Number for Rotor Blade Section. 0^ = 98** 
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Figure 24 — Concluded 
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Figure 25 Comparison of Theoretical Trailing Edge Noise Source Directivity 
Patterns in Measured and Retarded Coordinate Systems 
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Figure 27 One Third Octave Band Noise Spectra for NACA 0012 Data Reported 
by Schiinker, 0^ = 90*, a = 6®, rg = 2.81m 
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Figure 28 — One Third Octave Band Trailing Ed 
Data Reported by Brooks, =:90*', 
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Figure 29 



Normalized 1/3 Octave Band Trailing Edge Noise Spectra for Rotor Blade 
Section, 0- =98®, a = 7.6®, rfi = 3m 
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Figure 31 Normalized 1/3 Octave Band Trailing Edge Noise Spectra for NACA 0012 
Data reported by Brooks, =90^ a =5^ r^ = 1.22 m 
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a) THEORY USING FLAT PLATE SURFACE PRESSURE SPECTRUM 

Figure 33 Comparison of Measured and Theoretical 1/3 Octave Band Noise Spectra 
for Rotor Blade Section, 0c = 70®, a =7.6®, ro = 3 m 
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b) DIFFERENCE BETWEEN AIRFOIL AND FLAT PLATE SPECTRUM LEVELS 

Figure 35 Comparison of Airfoil and Flat Plate Normalized Surface Pressure 
Spectra, Bandwidth, 1 Hz 
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Figure 37 Comparison of Normaiized 1/3 Octave Band Spectrum Functions for First 
Principles Theory and the Scaling Law, ^9 sgo** 
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a) FLYOVER SPEED OF 72 m/sec, 6 q = 90 ° 



FREQUENCY, KHz 

b) FLYOVER SPEED OF 36 m/sec, 6»e = 90“ 


Figure 39 Comparison of Measured 1/3 Octave Band Helicopter Flyover Noise 
Spectrum and Predicted Trailing Edge Noise Spectrum 
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Figure 40 Diffraction Pattern of the Directional Microphone System for a 
Point Source of Sound at M = 0, R s 2.8m 




Figure 41 Gain and Diffraction Pattern Haif-width of Directionai Microphone 
System for a Point Source of Sound at M = 0 
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Figure 44 Locating the Directional Microphone Reflector to Account 
for the Apparent Source Position 












ACOUSTIC SOURCE 



Figure 45 Physics of Scattering 
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Figure 46 Variation of Acoustic Source Diffraction Pattern with Open Jet Mach 
Number, = 90** 







Figure 47 Gain and Apparent Diffraction Pattern Half-width of 
Directional Microphone System, M =0.5, = 90** 
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a) APPARENT SOURCE POSITION IN THE HORIZONTAL PLANE CONTAINING THE SHEAR 
LAYER NORMAL AND THE FLOW VECTOR 



b) APPARENT SOURCE POSITION IN THE PLANE NORMAL TO THE X-Y PLANE 
AND AT AN ANGLE 6 ^ WITH THE SHEAR LAYER 


Figure 48 Coordinates Defining Apparent Source Position 
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